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ABSTRACT 

We report near-infrared (primarily if-band) adaptive optics (AO) imaging 
with the Gemini-N and Subaru Telescopes, of a representative sample of 32 
nearby (z < 0.3) QSOs selected from the Palomar-Green (PG) Bright Quasar 
Survey (BQS), in order to investigate the properties of the host galaxies. Two- 
dimensional modeling and visual inspection of the images shows that ~36% of 
the hosts are ellipticals, ~39% contain a prominent disk component, and ~25% 
are of undetermined type. Thirty percent show obvious signs of disturbance. 
The mean M H (host) = -24.82 (2.1 Z£), with a range -23.5 to -26.5 (~0.63L£ 
to lOLJj). At < Lft, all hosts have a dominant disk component, while at >2L H; 
most are ellipticals. "Disturbed" hosts are found at all Mr (host), while "strongly 
disturbed" hosts appear to favor the more luminous hosts. Hosts with promi- 
nent disks have less luminous QSOs, while the most luminous QSOs are almost 
exclusively in ellipticals or in mergers (which presumably shortly will be ellipti- 
cals). At z < 0.13, where our sample is complete at 5-band, we find no clear 
correlation between M B (QSO) and M H (host). However, at z > 0.15, the more 
luminous QSOs (Mb < —24.7), and 4/5 of the radio-loud QSOs, have the most 
luminous if -band hosts (> 7 L^), most of which are ellipticals. Finally, we find 
a strong correlation between the "infrared-excess", Ltr/Lbb, of QSOs with host 
type and degree of disturbance. Disturbed and strongly disturbed hosts and 
hosts with dominant disks have L ir /L B b twice that of non-disturbed and ellip- 
tical hosts, respectively. QSOs with "disturbed" and "strongly-disturbed" hosts 
are also found to have morphologies and mid/far- infrared colors that are similar 
to what is found for "warm" ultraluminous infrared galaxies, providing further 
evidence for a possible evolutionary connection between both classes of objects. 

Subject headings: instrumentation: adaptive optics — galaxies:hosts — galax- 
ies:QSOs 
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1. Introduction 

It is now widely accepted that quasi-stellar objects (QSOs) are powered by accretion onto 
a massive black hole (MBH) within a galaxy. While MBHs are commonly found in galaxies, 
QSOs require the central black hole to be fed by a steady supply of matter. Characterization 
of the host galaxies of QSOs can provide important clues to understand both the formation 
mechanisms and the evolution of QSOs. 

Observations of low-z QSOs offer the obvious advantage of increased resolution and sen- 
sitivity for accurately measuring host properties, and relatively large samples are required to 
both adequately probe the full range of host properties and to explore possible relationships 
between host type and properties of the QSO nucleus. The pioneering optical surveys of low-z 
QSOs in the mid 1980's, using 4m-class telescopes (Hutchings, Crampton, & Campbell 1984; 
Smith et al. 1986), produced several common themes, as well as enigmatic problems that still 
remain today. QSOs were found typically to live in luminous (> L*) hosts, with a mixture 
of host types (spirals, barred-spirals, ellipticals), including a significant fraction of disturbed 
and/or strongly-interacting systems. There was also some hint of a trend in host type with 
host luminosity, with spirals and luminous ellipticals favored at the low and high end, re- 
spectively, of the range of QSO luminosities. However, it often proved difficult to distinguish 
reliably between different host types, often due to poor sensitivity or angular resolution or 
both; in addition, the mixture of objects being observed was abnormally weighted toward 
radio-loud and/or X-ray selected objects despite these objects being only a small fraction of 
the predominant radio-quiet population. Optical QSO surveys with HST/WFPC2 (Bahcall 
et al. 1997; Boyce, Disney, & Bleaken 1999; McLure et al. 1999; Hamilton et al. 2002; Dun- 
lop et al. 2003) offered dramatically improved resolution, but for the most part contained 
similar biases toward radio-loud objects, resulting in a preponderance of ellipticals. In the 
end, results from the HST surveys did not significantly modify earlier conclusions based on 
the ground-based surveys. 

More recent surveys of QSO hosts have favored the near-infrared, where the fraction 
of the observed light due to the underlying host galaxy is largest, and ground-based seeing 
can rival that from HST. Important large near-infrared QSO surveys include those from the 
ground with 2/4m-class telescopes (Dunlop et al. 1993; McLeod & Rieke 1994a,b; Taylor et al. 
1996; Percival et al. 2001) and, more recently, those exploiting adaptive optics (AO) (Marquez 
et al. 2001; Surace, Sanders, & Evans 2001), and, finally, those using HST/NICMOS (McLeod 
& McLeod 2001; Veilleux et al. 2005). Of most relevance to our own current work are the 
near-infrared observations by McLeod and co-workers who have imaged the host galaxies for 
a large volume-limited sample of optically-selected QSOs from the PGBQS similar to the 
sample and selection criteria we have used. Other near-infrared surveys of QSOs include 
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studies with nearly equal mixtures of radio-quiet and radio-loud QSOs (e.g., Dunlop et al. 
1993; Taylor et al. 1996; Marquez et al. 2001; Percival et al. 2001), and studies of "far- 
infrared-excess" QSOs (e.g. Surace, Sanders, & Evans 2001; Veilleux et al. 2005). All of 
the near-infrared surveys generally agree with previous optical surveys in finding that QSOs 
typically reside in luminous (i.e., > L*) hosts, and that spirals represent an increasing fraction 
of hosts at lower QSO luminosities. However, the near-infrared surveys still emphasize the 
overall diversity of host types, and they are still plagued by difficulties with the accurate 
determination of radial profiles due to a combination of problems including bad seeing, 
unstable PSF, and ultimately the faintness of the hosts at large radii. 

Finally, it has now become clear that optically-selected QSOs may be rivaled in their 
numbers by newly discovered populations of near- infrared selected (Cutri et al. 2002), and 
far-infrared selected (Low et al. 1988; Sanders et al. 1988) QSOs. It is extremely intriguing 
that a very high fraction (> 70%) of the former show evidence of tidal interactions (Hutchings 
et al. 2003), and nearly all of the latter are found to be major mergers of large gas-rich 
galaxies (e.g., Sanders & Mirabel 1996). These new findings raise the obvious questions of 
possible evolutionary connections between these different classes of objects, all of which seem 
to be capable of producing luminous QSOs. It thus became clear to us that it might prove 
extremely useful to revisit the properties of optically-selected QSO hosts and to search for 
trends in host type with other properties such as far-infrared luminosity and color, which 
might then prove useful in future studies designed to explore possible connections between 
the familiar optically-selected QSO population and the various types of infrared-selected 
QSOs. 

The availability of AO systems on 8-10 m class telescopes provides a potentially im- 
portant new capability for accurately measuring the distribution of light in the host galaxy 
of QSOs. In particular, our access to relatively large amounts of observing time using the 
new University of Hawaii near- infrared AO system, Hokupa'a, on Mauna Kea, allowed us 
to carefully design a sensitive and relatively uniform survey, primarily at i/-band, of a large 
representative sample of low-z QSOs randomly selected from the optical (i?-band) PGBQS. 
The sample was made large enough to include a representative sub-sample of the minority 
population of radio-loud objects, as well as to allow for adequate sampling of a wide range of 
QSO luminosities, but was otherwise unbiased. An additional bonus for our new observations 
was the existence by 2003 of relatively complete and accurate spectral-energy-distributions 
(SEDs) from radio-to- Xray wavelengths (Sanders et al. 1989; Elvis et al. 1994), including 
most importantly, new far-infrared fluxes for most of our targets, which were recently deter- 
mined from reprocessed IRAS and ISO data (Haas et al. 2000, 2003). 

In this paper we present the first near-infrared imaging study of a large number of 
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QSOs (32 objects) made with AO systems on the Gemini 8.2m and Subaru 8.1m telescopes 
on Mauna Kea. We first describe, in §2 how our sample was selected. The observations 
and data reduction techniques are detailed in §3. The QSO host galaxies are characterized 
through 2-D fits and visual analysis of the images in §4. Section 4.4 presents our adopted 
classification scheme, including descriptions of the individual objects and comparisons to the 
host types derived from previous imaging studies, and summarizes the host fit parameters. 
Section 5 discusses how host type depends on other properties of the QSO (e.g. degree 
of disturbance, radio loudness, the near-infrared absolute luminosity of the host, and the 
total far- infrared luminosity), and compares our results with previous studies. Our main 
conclusions are summarized in §6. 

As background to our current work, good brief reviews of previous problems and sub- 
tleties associated with interpreting the host types of low-z QSOs can be found in McLeod 
& Rieke (1994a) and Marquez et al. (2001). In addition, a more thorough review of the 
wide range of QSO host studies - their consensus and differences - including theories as to 
the origin and evolution of QSOs and their possible relation to other types of extragalactic 
objects, can be found in the proceedings of three recent workshops specifically devoted to 
the study of QSO hosts (Clements & Perez-Fournon 1997; Marquez et al. 2001; Barthel & 
Sanders 2006). 



2. Sample selection 

The performance of AO systems is a strong function of the photon flux seen by the 
wavefront sensor (WFS). Although the Hokupa'a (the Institute for Astronomy AO system 
(Graves et al. 1998) in use on the Gemini North Telescope from mid-2000 to mid-2002) and 
Subaru AO (Takami et al. 2004) systems can "lock" onto objects as faint as my ~ 18 under 
good seeing, only very partial wavefront correction can be achieved for such faint objects. 
To ensure reasonably good image quality (FWHM< 0'/2), our sample was limited to bright 
QSOs (approximately m v < 16). This observational constraint would result in a strong 
bias towards intrinsically highly luminous QSOs if the sample were not redshift-limited. 
Moreover, observation of higher redshift QSOs is made more difficult for the purpose of this 
work by the decrease in spatial resolution, and especially surface brightness, with increasing 
distance. Our sample is hence limited to z ~ 0.3. 

The ideal sample for our current observations proved to be the Palomar-Green Bright 
Quasar Sample (BQS: Schmidt & Green 1983). Although there are more recent and larger 
QSO surveys, the BQS still contains the majority of the nearest and brightest optically- 
selected, UV-excess QSOs for study from the northern sky, and its magnitude limit was 
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close to our desired brightness limit for good AO performance. There are also extensive 
multi-wavelength observations in the literature for the BQS sample, which potentially allows 
us to further explore the relationship of the QSO host properties to other properties of the 
QSO. 

The BQS contains 92 QSOs with M B (Vega) < -23 (in the H a = 50, q Q = cosmology 
used in Schmidt & Green 1983), in the 10,714 deg 2 surveyed region of the sky at |6| > 30° 
and 5 > —14°. Forty-seven of these QSOs met our final redshift and brightness selection 
criteria (z < 0.3, By cga , < 16.3), of which we were eventually able to observe 58% (28/47) 
during our two year observing campaign on Mauna Kea. These 28 QSOs from the BQS were 
essentially selected at random from the complete sub-sample of 47, depending only of the 
RA range observable during each observing run. 

Four additional objects were eventually added to the list of 28, bringing the total number 
of observed targets to 32. Two objects from the BQS with M B > -23 (PG1119+120 and PG 
1126-041) were successfully observed during one of the earliest observing runs, while we were 
still defining the scope of the overall AO campaign. Rather than discard the data, we decided 
to include these objects here for two reasons. First they allow us to further explore the low 
luminosity end of the QSO sample, and second, we noted early on that the uncertainties 
in the S-band magnitudes listed in the original BQS are large enough such that these two 
objects might have properly been in our sample. Two additional objects were added to our 
observing list in the latter half of our overall campaign, when it was realized that we had only 
two "bona-fide" radio-loud QSOs (RLQs) in our list of observed objects. One of these RLQs, 
4C31.63 [z = 0.295), was at the time it was added, the highest redshift object observed, but 
still within our planned redshift limit for the survey, and its _B-band magnitude was above 
the average overall threshold for the BQS. The second RLQ, PG2251+113 (z = 0.326), was 
above our stated redshift cutoff, but satisfied our 5-band limit and was the only RLQ from 
the BQS that could be observed at the RA range available for observations during 2001, 
September. The addition of these two RLQs allowed us to bring the total fraction of RLQs 
in our observed sample (5/32 ps 16%) up to the ratios more typically found (i.e., ~15%) in 
large QSO surveys. 

The complete list of our 32 observed targets is given in Table 1 (Vega magnitudes are 
used in this table and throughout this paper). The listed _B-band magnitudes have been 
updated where possible with newer more accurate values than originally given in Schmidt & 
Green (1983). The absolute magnitudes are also updated to an H = 75, Qm — 0.3, Q\ = 0.7, 
cosmology: the M B < —23 QSO criteria used by Schmidt & Green (1983) now corresponds 
to M B < —22.1. We note that the updated S-band magnitudes resulted in three additional 
objects, PG0838+770, PG0844+349 and PG1229+204, falling below the absolute magnitude 
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cutoff, bringing the total number of lower luminosity objects to five. 

3. Observations and data reduction 
3.1. Instruments 

Both the Hokupa'a AO system, built by the University of Hawaii and used on the 
Gemini North telescope, and the Subaru Telescope AO system, were used for this project. 
For both systems, a 36-element curvature sensor is used with a 36-element bimorph mirror. 
Curvature AO systems are traditionally very efficient on faint targets, partly thanks to the 
use of photon-counting devices (Avalanche Photo Diodes) with no readout noise. This allows 
these systems to guide on my = 18 sources and to routinely provide good correction (better 
than 0'/2) on m v = 16 point sources, which makes them well suited for QSO host galaxy 
imaging. The AO-corrected beam was fed to QUIRC for Gemini and to IRCS for Subaru. 
Both cameras use 1024 x 1024 detectors for near-IR imaging. 

3.2. Choice of filter and target 

The contrast between the QSO and a typical host galaxy is better (lower) at about 1 
(Km (in the rest frame) than either in the visible or the longer infrared wavelengths. This is 
largely due to a local minimum near 1 /im in the rest-frame radio-to- Xray spectral energy 
distributions (SEDs) of nearly all QSOs, both RLQs and RQQs, except for the strongly 
beamed Blazars — e.g., Fig. 2 of Sanders et al. (1989), and also Fig. 1 of McLeod & Rieke 
(1995). Moreover, the performance of AO systems is optimal in the near infrared: at shorter 
wavelengths, the number of actuators is too small to sample the wavefront properly. At 
longer wavelengths, the high thermal background precludes the observation of low-surface- 
brightness sources, and the diffraction limit of the telescope restricts the achievable angular 
resolution. Observations in the near-IR are therefore well-suited for QSO host studies from 
the ground. 

The seeing was constantly estimated during the observations by image quality monitor- 
ing and monitoring of the control signals of the AO system. Most of the observations were 
made using the H filter, which provides the optimal resolution under good seeing conditions. 
When the seeing was bad (corrected FWHM in H band larger than 0V2), the K' filter was 
preferred to the H filter, and usually allowed high resolution imaging (FWHM< 0'.'2) even 
under worse than average seeing, but at the expense of a stronger thermal background. A 
few observations were done in the J band on exceptionally good nights. Fainter QSOs were 
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preferentially observed under good seeing conditions, to avoid running out of bright targets 
when the seeing was average or worse than average. These precautions resulted in a surpris- 
ingly homogeneous set of images (angular resolution about 0'.'15) considering the wide range 
of seeing conditions and QSO magnitudes of the sample. 

Table 2 gives the complete summary of our observing log for all of the data that are 
presented in this paper, including filter, exposure times, seeing (FWHM), and telescope used. 
Twenty-six objects were eventually observed at H, and 8 of these 26 were also observed at 
either J or K', or both. Four objects were observed at K' only, one at J only, and one 
at both J and K'. Twenty-eight objects were observed at Gemini, and four objects were 
observed at the end of our campaign (Feb, 2003) using Subaru. 

3.3. Reference PSF observations 

For each target, a list of about six PSF candidates, all within a few arcminutes of the 
quasar, was prepared prior to observations. Each of those was briefly observed (without 
closing the AO loop) to test if it could, with the help of the WFS neutral density filter set, 
reproduce the WFS photon counts generated by the quasar. This allowed us to match the 
photon counts of the quasar to within 10% for most of the objects in our sample. The AO 
loop parameters (membrane stroke and loop gain) were not modified between observations 
of the target and the PSF. 

For most of the targets, the PSF total integration time was much smaller than the QSO 
total integration time. Because of the variability of the PSF presented in §3.4, the residuals 
of the PSF subtraction are limited by effects that do not average well with time, and there is 
little gain in observing the PSF star for longer than a few minutes. However, long exposure 
times are required to record the faint structure of the host galaxy that are outside the central 
part of the PSF, where PSF subtraction is not essential (where the PSF wings are fainter 
than the host galaxy). When possible, the PSF star was observed with the same frame 
integration time as the QSO to use the PSF frames for sky extraction. 

3.4. PSF stability — lessons learned 

In this section, we characterize the PSFs delivered by Gemini+Hokupa'a and the Subaru 
AO system and investigate their variability. This preliminary work is essential for attempting 
the characterization of the QSO hosts in our data. 

The PSF delivered by an AO system is a function of : 
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• 1. Guide star V magnitude and color of the object imaged. The brightness 
of the star, as seen by the WFS, strongly affects the performance of the AO system. 
Due to the wavelength-dependence of the seeing (and the optics of the camera), PSFs 
are also expected to be color- dependant. 

• 2. Atmospheric seeing. Under poor seeing conditions, the FWHM of the corrected 
PSF is larger. In addition, the brightness of the large halo surrounding the PSF core 
increases with the turbulence. 

• 3. Instrumental variables. Wind-powered telescope shake, change in the optical 
path, and the primary and secondary mirror wavefront all modify the PSF. 

• 4. AO loop parameters. Loop gain and membrane stroke (for curvature systems) 
set the speed of the AO loop and the optical gain of the WFS respectively, and therefore 
affect the PSF. 

Despite efforts to match the brightness (as seen by the WFS) of the reference PSF and 
the QSO, small discrepancies in the two flux levels are unavoidable. More importantly, we 
found that the variations in seeing always limit the accuracy of the PSF subtraction when 
the flux of the PSF and the QSO are matched to better than 10%. Figure 1 shows how 
the radial profile of the PSF can change rapidly (timescale of about a minute) due to seeing 
variations. The far wings of the PSF (1" to 2" in Fig. 1) show how much light is scattered 
by small-scale turbulence (high spatial frequencies of the seeing), which is not corrected by 
the AO system. In this figure, it can be seen that these wings are much higher for PSF1 
than for PSF2. Consequently, because the turbulence was stronger during the acquisition 
of PSF1, the AO correction is poorer and, inward of about (X'2 (dotted vertical line), the 
surface brightness of the PSF is lower. The "crossing point," for which the loss of surface 
brightness due to the poor AO correction exactly equals the gain of surface brightness due to 
the increased scattering by small-scale turbulence, is set by the number of actuators of the 
AO system. The surface brightness at this crossing point is insensitive to seeing variations. 
Figure 1 also shows how the brightness of the PSF affects its radial profile: for brighter PSFs, 
the AO correction is better and flux is taken from an "intermediate" annulus between O'/l 
and (X'3 and put into the central diffraction core of the PSF. The far wings of the PSF are 
left unchanged because of the limited number of actuators of the AO system. During real 
observation, a combination of PSF variations due to changing seeing and different source 
brightness will affect the radial profile. 

Although the observer has no control over the seeing, it is possible to reduce PSF vari- 
ations by avoiding any change in the instrument during observations. For example, the AO 
loop parameters should be kept constant. During this observing program, it was discovered 
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that the PSF stability could be significantly improved by turning off the Cassegrain rotator 
of the Telescope. At the Gemini and Subaru telescopes, which have alt-azimuth mounts, the 
Cassegrain units carrying the AO systems are usually rotated to preserve the field orienta- 
tion on the detector. The primary and secondary mirror create static wavefront aberrations. 
If these static aberrations are rotating in front of the AO system, which has preferential 
directions due to the geometry of the actuators and also introduces some higher order aber- 
rations, the resulting rotational shear between the incoming wavefront and the aberrations 
of the AO system changes the structure of the PSF. This very complex PSF variation is 
difficult to predict since the individual sources of aberrations are poorly known. Figure 2 
illustrates this effect through the comparison of two differenced PSFs, one for which the 
Cassegrain rotator angle was kept unchanged between the two exposures, and the other for 
which a 15° rotation occurred between the two exposures. Because some of the aberrations 
come from the primary and secondary mirror and some come from the AO system and the 
camera, a PSF obtained after rotation of the Cassegrain rotator cannot be reproduced by a 
rotation of a PSF obtained at a different Cassegrain rotator angle. 



3.5. Decoupling of the PSF and the host galaxy 

We have seen above that turning off the Cassegrain rotator makes the PSF more stable 
and therefore, improves the accuracy of the PSF subtraction. There is a second advantage 
to this technique: the PSF and the image of the host galaxy can be decoupled through 
differential rotation. 

When attempting to detect faint structures very close to the QSO, we are usually very 
sensitive to small variations in the PSF structure, such as "flares" of the PSF. These "flares" 
are more frequent when observing faint targets such as QSOs, because the limited number of 
incoming photons doesn't allow fast recovery from events such as a wind-induced telescope 
shake. It is possible to encounter such a flare on the object and not on the reference PSF, 
leading to the false detection of a structure in the host galaxy. Careful visual examination of 
individual frames can help detecting these artifacts but is not fully reliable. When observing 
with the Cassegrain rotator off, the PSF has a constant orientation on the detector while the 
object rotates, which unambiguously decouples the object from the PSF variations. For a 
"flare" to be misidentified as a structure in the host galaxy, it would need to rotate in the PSF 
at a speed reproducing the field rotation. Because we observed most objects continuously for 
more than an hour, the total rotation during the observation is usually sufficient to efficiently 
decouple the object from the PSF. 

Because of the improved PSF stability when observing without the field derotator, most 
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of our observations were done with the Cassegrain rotator turned off. The geographical 
coordinates of the telescope, the UT time of the exposure, the coordinates of the target 
and the angle of the Cassegrain rotator are required, for each exposure, to recover the field 
orientation. 

This observing mode requires the exposure time to be short enough to avoid rotational 
smearing of the image. The limit we set on the rotational smearing corresponds to half a 
full width at half maximum (FWHM) at the edge of the detector when the target is at the 
center. This limit allows exposure times of 20 s to 30 s on most of the sky. 

For some objects, detector saturation was reached in a few seconds. Because of the 
relatively long readout time of the camera (10 s) on Gemini/Hokupa'a, we chose to obtain 
a series of short unsaturated exposures and a series of longer exposures which were essential 
to reach a good flux sensitivity. We were careful to obtain the longer exposures when the 
field rotation was slow. 



3.6. Host galaxy characterization with AO images 

AO PSFs are time-variable and have relatively strong and extended wings compared 
with Hubble Space Telescope (HST) PSFs, but the image of the underlying host is made 
sharp, which allows accurate retrieval of the radial profile in regions where the wings of the 
central PSF are not significantly brighter than the host galaxy. Depending on the objects of 
our sample, we find that it is possible to accurately recover the host radial profile outward of 
about 0'. / 2 (for hosts which account for a large fraction of the total flux) to 2" (for hosts which 
are much fainter than the central point source) with classical PSF subtraction techniques. 
This corresponds to about 0.5 kpc to 2 kpc for most of the objects in the sample. Thanks 
to the large collecting area of 8 m class telescopes, the host radial profile can usually be 
recovered to up to 5" to 10". Beyond this angular separation, the error due to imperfect sky 
subtraction becomes larger than the host: due to the small field of view of the images, the 
sky level cannot be accurately decoupled from the faint outer part of the host galaxy. 

The good angular resolution (O'/l to 0'.'2, which is comparable to or better than that of 
HST at the same wavelength) and depth of QSO images obtained on a 8 m class telescope 
equipped with an AO system makes it possible to image small and faint structures in the 
host. This is very important for this project, which relies on the ability to characterize the 
morphology of the host beyond the disk vs. elliptical galaxy classification. 
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3.7. Image coaddition and PSF synthesis 

For each target, the frames were processed and coadded by a data reduction pipeline 
written in C. Flat field and bad pixel maps were created from a set of dome flat exposures 
with increasing exposure time. The sky frame was obtained from the dithered QSO and 
nearby PSF exposures. 

The persistence effect of the detector was corrected for by ignoring the pixels that have 
been recently strongly exposed to light. An empirical formula was used to select such pixels 
from the pixel value and the time of the exposure of the 10 preceding frames. 

The first steps described above (flat fielding, bad pixel identification and sky subtraction) 
were applied to both the PSF and the QSO frames. For each QSO, a coadded image was 
generated by derotating each QSO frame before coaddition. As shown in Fig. 3, because of 
the combination of dithers, rotations and bad pixels masking (or pixels ignored due to the 
persistence of the detector), the total exposure time varies across the final image. 

A reference PSF frame is obtained through coaddition of individual PSF frames, but 
without field derotation. To avoid increasing the noise of the PSF-subtracted image, the far 
wings of the PSF were replaced by synthetic PSF wings generated from a fit of the wings 
radial profile. The PSF image matching the coadded QSO image was obtained by coadding 
multiple copies of the reference PSF frame with the proper rotation angles and relative 
exposure times. 

4. Analysis and Results 
4.1. Images 

Images are shown in figure ?? for each object of our sample. The flux in the raw images 
(left) is dominated by the bright central nucleus, although the host galaxy is usually visible 
without PSF subtraction, thanks to the good contrast and high angular resolution delivered 
by AO. The PSF-subtracted images (middle) reveal more host galaxy details, especially in 
the central part of the image. 

The quality of the PSF subtraction is highly variable through the sample, and depends 
upon many factors (see §3.4). In some cases, the reference PSF used for subtraction is less 
peaked (and, therefore, has a comparatively brighter "halo" component) than the actual 
PSF during observation. If the host galaxy vs central nucleus luminosity ratio is low, the 
residual of the PSF subtraction is then negative in an annulus surrounding the positive resid- 
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ual at the center of the image (see §4.3 for a description of the 2-D analysis algorithm used 
to subtract the PSF): this effect can be seen on PG0026+129, PG0838+770, PG1411+442. 
PSF subtraction errors due to poor stability of the AO control also occur when the source is 
faint, the atmospheric turbulence is strong and/or clouds vary the source brightness during 
the observation. These effects are responsible for the poor PSF subtraction of PG0050+124, 
PG1116+215 and PG1302-102. During the observation of PG1617+175, an avalanche pho- 
todiode of the AO wavefront sensor burned out, suddenly creating a bright extension in the 
PSF delivered by the system, as can be seen in the PSF subtracted image. 

Many of the images also show artefacts due to sky subtraction residuals. These fea- 
tures are usually faint and beyond 5" from the QSO. They are most visible in the images 
of PG0804+761 (curved feature 4" North- West of the QSO and compact artefacts 10" due 
East, North, West and South of the QSO), PG0838+770 (about 5" to 7" South of the 
QSO), PG0844+349 (straight edge and extended feature 6" East of the QSO), PG0923+201 
(straight feature running North-South about 5" East of the QSO), PG1001+054 (faint ex- 
tended fringes 5" to 10" South of the QSO), PG1119+120 (faint extended fringes mostly in 
10" South- West of the QSO), PG1307+085 (4 faint features in a 10" square pattern centered 
on the QSO), PG1626+554 (wide 5"-radius arc centered in the QSO, mostly visible in the 
northern half of the image) and PG2214+139 (faint broad features in a 8" square pattern 
centered on the QSO). 



4.2. Photometry 

The flux within a 5"5 diameter circular aperture centered on the QSO was measured 
for each object, with an accuracy of about ±5% at the 1 a confidence level. This aperture 
size allows easy comparison with previous measurements by Neugebauer et al. (1979) and 
Neugebauer et al. (1987). For the objects observed in non-photometric conditions, only a 
lower limit for this flux could be determined. 

Table 3 compares our photometry with other measurements from the literature and 
from the 2MASS archive. No correction was applied to account for differences between the 
2MASS K s and Mauna Kea K' photometric systems, as these two filters are very similar 
(Tokunaga et al. 2002). 

The values shown in Table 3 are consistent with previous findings that show that many 
of the QSOs observed in this sample appear to show some variability in the near-infrared. 
The amplitude of variability shown in this table is in good agreement with Neugebauer et al. 
(1989): both sets of data identify PG1001+054, PG1426+015, PG1613+658, PG1617+175 
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and PG2 130+099 as the most variable QSOs of our sample. Since the observations made 
under non-photometric conditions need to be calibrated with other photometric measure- 
ments of the same object, the variability of QSOs can cause errors in the estimate of the 
host galaxy luminosity. In order to minimize such errors, the time difference between our ob- 
servations and the photometric measurement should be kept small. 2MASS measurements, 
rather than older measurements from Neugebauer et al. (1979) or Neugebauer et al. (1987), 
were therefore used for photometric calibration of objects for which no reliable photometry 
could be derived with our observations. In one case (PG1116+215), the lower limit obtained 
with our data is higher than the 2MASS measurement: this lower limit was then adopted for 
photometric calibration. We note that PG0050+124 was observed under very poor condi- 
tions: poor seeing and large absorption by clouds resulting in highly variable adaptive optics 
performance. We did not attempt to perform radial profile or x 2 analysis for this object, 
which explains why no photometry was adopted in Table 3. 



4.3. Two-dimensional x 2 analysis of the hosts 

The 2-D analysis relies on a least-square fitting of the image according the the following 
equation: 

2 (Image — (Point Source + Host) ® PSF) 2 ^ 

pixels 

Pixels too close to (distance less than 1" to 2" depending on the QSOs), or too far from 
(more than 8" to 10") the center were excluded from the fit because PSF calibration errors 
dominate the host's signal in the inner region of the image, and sky calibration is unreliable 
in the outer parts of the image (due to the small field of view of the detector). Rejecting 
the central part of the image in the fitting algorithm could allow solutions where the total 
flux of the image is not preserved: the PSF contribution could be increased to a point where 
its total flux exceeds the total flux in the raw image to fit the host galaxy beyond the inner 
edge of the fitting region. To avoid this effect, the total flux of the PSF component was not 
allowed to exceed the total flux in the raw image. 

Bright companions, when present, were also excluded from the fit. The value of a was 
directly measured on the image. This 2-D fit can have up to 9 free parameters: 1 for the 
intensity of the central source, 4 for the bulge component (modeled in this work as a single 
de-Vaucouleurs r 1//4 profile) and 4 for the disk component (modeled as an exponential profile). 
A simulated annealing minimization routine was used to find the absolute minimum. 

The x 2 minimization routine was first run with the full 9 free parameters to yield 
X^ xp+r i /4 , the minimum value of x 2 f° r a exponential + r 1 / 4 de Vaucouleurs profile host 
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model. Then, ( r1 ^ 4 nost on ly) an d Xe xp ( ex P nos t onr y) were computed, with 5 free 

parameters in each case. The resulting reduced x 2 values for each of these three cases are 
listed in Table 4 for all of our targets. It is worth noting here that the noise component a 
in equation 1 is often poorly characterized, due to the difficulty of estimating the temporal 
variations of the PSF wings. Thus the reduced x 2 values listed in Table 4 cannot be reliably 
compared between targets, but different fits obtained from the same image can still be 
compared to classify the host galaxies. 

For each object, the 1-D radial profiles extracted from the 2-D image and the corre- 
sponding 2-D fit models (point nucleus + host model convolved by the PSF) are shown in 
Fig. 5 to ??. The 1-D profile of the residual error, which is also shown in the same figures 
for the 3 fits, tends to become large at small angular distance from the center because of 
errors in the PSF subtraction. The domain within which the fit was performed is shown by 
2 vertical red lines. 

Due to the limited field of view of the instrument (20"), estimations of sky background 
level and the PSF are prone to errors beyond ^8" radius. This explains the peculiar behavior 
of some profiles beyond ^8" (for example, the abrupt drop at 10" in the PSF profile for 
PG1440+356), and justifies our choice to exclude the outer part of the profile for host 
galaxy characterization. 

4.4. Adopted host classifications 

4-4- 1- Classification scheme 

In deciding which of the three x 2 - m odel fits listed in Table 4 best characterizes the host 
galaxy, the following criteria were used. Unless the 2-D fit strongly favors a 2-components 
bulge+disk host (i.e. if X 2 i/4 +cxp < m ^ n (Xcxp' X^i/*) — 0-1 an d no more than 70% of the total 
host flux is allocated to either the bulge or disk in this bulge+disk fit), the smallest of the 
X 2 i/4 and xl X p is adopted as the result of the 2-D fit. The adopted fits for each object and 
each filter used are listed in columns 2-4 of table 5. 

However, in a number of cases, the difference in the x 2 values for the exponential vs. the 
r 1//4 de Vaucouleurs profile were small enough to render the choice of a particular host model 
somewhat meaningless. In the cases where 0.02 < \x 2 i/4~xl xp \ < 0-05, we consider the result 
of the fit to to favor one model, but with very low confidence (noted in parentheses in table 
5). And in those cases where \x 2 i/4 — Xc xp l ^ 0.02, or in the case of bad PSF subtraction, we 
consider the fit to be inconclusive (a question mark in table 5). 
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Since the overall QSO host classification scheme used in this work takes into account both 
the results of the 2-D x 2 -model analysis and morphological information visually extracted 
from the images, table 5 also lists information obtained from visual inspection of the host 
images, including the presence or absence of "Bars," (spiral) "Arms," and the overall degree 
of "Dist(urbance)," respectively. Additional visual information is given in the Notes. 

Each host was finally classified using a combination of information which included the 
best-fit x 2 - m odel "type" and the visual appearance. The final "Adopted Type" listed in col. 
8 of table 5 was either "D p " = Disk present, "E" = Elliptical, or "B+D" = Bulge+Disk, 
except for those objects with insufficient and/or ambiguous information from the x 2 ~ m odel 
fits, or objects with "strongly disturbed" disks, where a "?" was listed. We note here that 
since the 2-D fitting analysis does not include the central part of the image (due to PSF 
variability), the type "Disk present" does not exclude the presence of a bulge (which would be 
hidden in the central part of the image), and the "Disk present" /'"Bulge+Disk" distinction 
is therefore unclear. "Elliptical" and "Bulge" (in the Bulge+Disk type) refer to a r 1 / 4 for 
all 2-D fits performed in this work. Strongly disturbed hosts were not classified according to 
host type, since such a classification would be difficult and unreliable. 

While the degree of disturbance is estimated purely from visual inspection of the images, 
the rules adopted for the host "adopted type" classification were as follows: 

• If a bar or spiral arms are obvious in the image, the host is classified as a "Disk present" , 
or a "Bulge+Disk" (the decision between "Disk present" and "Bulge+Disk" is made 
using the 2-D fitting analysis results). 

• If no spiral structure or bar is visible in the image, the host classification is given by the 
2-D x 2 -model fitting analysis. If the host has been imaged in several spectral bands, 
and the corresponding 2-D x 2 -model fitting analysis do not agree, the fitting result 
obtained for the spectral band with the best PSF subtraction is adopted. 

It is worth restating here that we have chosen to give priority to the features visible 
in the images rather than the results of the 2-D x 2 -model fits, primarily because the 2-D 
X 2 -model fits are prone to errors. It can be seen from the 1-D profiles in Fig. 5 to ?? that 
for many objects, the 2-D fit does not produce a reliable results - we found that changing 
the fit domain often changes the result of the fit. Thus, it was reasoned that, for example, 
in several objects where a bar or spiral arms were clearly visible in the host images, the host 
should not be classified as an elliptical, even if the 2-D fit suggests so. It also should be 
noted that not all objects gave ambiguous results. Some objects are quite easy to classify, as 
the x 2 -model fit results are clear and agree with the images. But overall, the images show 
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that the x 2 " m odel fits often give misleading results, and in ~19% of the objects (6/32) the 
hosts are so disturbed that a x 2 -model fit is of little relevance. 

4-4-2. Notes on individual objects 

Before presenting the final fit parameters corresponding to the host "adopted type" 
listed in column 8 of Table 5, it is worth comparing our x 2 - m odel fit results and observed 
visual features with results from previous published studies of these objects. Each object is 
discussed separately in the paragraphs below, including comparisons with previous analyses 
by different authors and their published "host types," and including remarks about the 
various difficulties encountered in our own analysis, which often led to ambiguities in the 
final choice of host type. 

• PG 0026+129 

The image shows an obvious host, slightly elongated east- west, with no apparent struc- 
ture. The radial profile analysis fails to identify clearly the galaxy type, but gives host 
magnitudes ranging from m H = 13.43 (bulge+disk fit) to m H = 13.69 (disk fit), and 
nucleus magnitudes ranging from m H = 13.02 to m H = 13.05. Although our image 
does not allow classification of this host, the ifST/NICMOS observations of McLeod 
& McLeod (2001) suggest that it has an r 1//4 -law profile. A non-stellar apparent com- 
panion is visible 4.8" from the central nucleus (m H = 18.5, PA = 22°). 

• PG 0050+124 

This object is also known as I Zw 1 and is often classified as a Seyfert galaxy, partic- 
ularly since the host shows clear spiral arms (Bothun et al. 1984; Surace et al. 1998; 
Zheng et al. 1999; Canalizo & Stockton 2001). In our observations, the spiral arms are 
visible in the host galaxy across the field of the detector (20"). Although the observ- 
ing conditions for this object were poor (making it impossible to perform an efficient 
PSF subtraction or photometry), the morphology is unambiguously that of a spiral 
host. Canalizo & Stockton (2001) have found spectroscopic evidence for continuing 
star formation in the disk as well as morphological evidence for weak interaction with 
the small companion to the west. Scoville et al. (2003) find a mass in molecular gas 
of ~ 10 10 M Q , and recent observations with the BIMA interferometer (Staguhn et al. 
2004) resolve individual giant molecular clouds in a ring structure from about 1.2 to 
3.5 kpc from the nucleus. 

• PG 0157+001 

The host galaxy of PG 0157+001 (= MRK 1014) is large and bright, and it shows 2 
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prominent spiral-like features that extend from < 1 kpc to at least 28 kpc (limited 
by the edge of the detector, 10" from the nucleus), as well as a bright bulge. This 
galaxy has been the subject of numerous ground-based and HST imaging studies (e.g., 
MacKenty & Stockton 1984; Stockton & MacKenty 1987; Hutchings & Crampton 1990; 
Surace et al. 1998; McLeod et al. 1999; Canalizo & Stockton 2000). The spiral features 
are clearly not conventional spiral arms (see Fig. 2 of Canalizo & Stockton 2000); 
instead, as many of the previous studies have emphasized, they are almost certainly 
tidal tails generated in a past major merger. These tails provide strong evidence 
that both of the merging galaxies themselves had cold stellar disks. The host galaxy 
is highly asymmetric, with the northern tail brighter than the southern tail. Several 
bright knots can be seen in the northern tail, which are also clearly seen in HST images 
(e.g., McLeod et al. 1999). A small, bright companion is visible to the south-west at a 
projected distance of about 25 kpc. The host is well fitted as a mn = 13.30 bulge+disk 
galaxy with a mn = 13.41 central nucleus, which is compatible with the structure 
seen in the image. Canalizo & Stockton (2000) confirm the presence of significant 
populations of young stars, particularly along the leading edge of the northern tail. 

• PG 0804+761 

This host is best described as a "bar" about 5 kpc in radius. An extended (kpc size) 
faint companion is visible to the East of the QSO, at a projected distance of 4 kpc 
from the QSO. The curved structure visible in the image about 8" to the North- West 
of the QSO is a artifact from the camera. No reliable 2-D fit was found for this galaxy 
in K' ', and the fit in H suggests an elliptical galaxy. However, the bar indicates that 
this galaxy contains a prominent disk. The only previous published imaging study of 
this object is that of McLeod & Rieke (1994a), who found a profile consistent with 
a disk galaxy from their ground-based if-band imaging. Scoville et al. (2003) detect 
molecular gas mass of about 3 x 10 9 M in this galaxy. 

• PG 0838+770 

The central part of the host of PG0838+770 is a 9 kpc radius bar. A noticeable twist 
of this bar is seen in our image, and a 20 kpc long arm is starting from the East end 
of the bar. This arm seems to connect both ends of the bar, thus forming the northern 
half of a ring immediately surrounding the bar. Visible images (Surace, Sanders, & 
Evans 2001) show both the bar and two spiral arms outside the ring. An extended 
(about 2.5 kpc diameter) apparent companion is also visible to the northwest at a 
projected distance of 12 kpc. Although the host is best fitted by a r 1 / 4 law, the bar 
and spiral arms show that a substantial disk is present. 
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• PG 0844+349 

Two symmetric spiral arms are clearly visible, extending to about 9 kpc radius in our 
image. Hutchings & Crampton (1990) describe the host as a "barred spiral, with some 
complex outer structure." Our 2-D fit, which suggests a dominant bulge, is uncertain 
because of poor PSF subtraction. 

• PG 0923+201 

Although some previous studies have listed this object as interacting (McLeod & Rieke 
1994b) or possibly interacting (Hutchings, Janson, & Neff 1989), apparently because 
of the bright galaxy projected ~ 10" to the south, the host is unambiguously a very 
symmetric elliptical. It has an apparent extended companion at a projected distance 
of 7 kpc (272) almost due north. This host was also identified as an elliptical from 
HST/WFPC2 observations (McLeod et al. 1999). The 2-D fit shows that the host 
galaxy (mj = 15.34, m H = 14.73) has a low ellipticity (about 0.1), and the central 
nucleus is about 1 magnitude brighter than the galaxy. 

• PG 0953+414 

This host was identified as a spiral from ground-based imaging by Hutchings, Janson, 
& Neff (1989); however, all subsequent imaging studies, including HST/WFPC2 ob- 
servations of Bahcall et al. (1997) and McLeod et al. (1999), have found the host too 
faint to classify reliably. This is also true for our data. The brightest part of the host is 
visible as a faint extension to the south of the QSO in our image. McLeod et al. (1999) 
find marginal evidence favoring an elliptical profile, but they also claim a possible tidal 
arm to the southwest, which would indicate the presence of a disk. 

• PG 1001+054 

Our image shows a highly elongated host (long axis radius of 4 kpc) which is probably a 
bar but could also be a small inclined disk. The 2-D fit suggests an elliptical galaxy, but 
we choose to classify the host as a "disk present" from visual inspection of the image. 
Surace, Sanders, & Evans (2001) see only a small, uniform host from ground-based 
optical and near-IR imaging. 

• PG 1004+130 

The quasar, also known as 4C 13.41 and PKS 1004+13, is a powerful steep-spectrum 
radio source. The host galaxy is large, bright and smooth, and is shown to be an 
elliptical by the 2-D fit. Bahcall et al. (1997) also classify the host as a bright elliptical 
from HST/WFPC2 imaging. 

• PG 1116+215 

The PSF subtraction is poor for this object. Nevertheless, the host appears to be 
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asymmetric, with several companions: a compact companion at a projected distance 
of about 4 kpc to the west, a faint extended (2 kpc diameter) companion at 5.5 kpc 
projected distance to the south, a compact companion (star?) at 12.7 kpc to the south- 
southwest, and an elongated bright galaxy at 33 kpc to the West (outside the field of 
view of the image shown in Fig. 5). The 2-D fit to the host suggests that it is a disk, 
although there is considerable uncertainty in the determination. Bahcall et al. (1997) 
list the host as "probably an elliptical." 

• PG 1119+120 

The host is bright and smooth, with a north-south elongation. An elongated companion 
is visible 2.6 kpc (2'.'7) to the north, projected on the host galaxy. Optical images 
(Surace, Sanders, & Evans 2001) show 2 faint extended rings on the east and west 
sides of the brightest part of the host. These features are too faint to be detected in 
our Near-IR images. Our 2-D fitting gives a bulge+disk decomposition and shows that 
the host galaxy is almost a magnitude brighter than the central nucleus. 

• PG 1126-041 

Our image shows a bright and large host with spiral arms extending beyond 10 kpc. 
The 2-D fit suggests the presence of a bulge of luminosity comparable to the disk 
component. A bright apparent companion, at a projected distance of 6.6 kpc (5'.'7), is 
clearly resolved and shows a small elongation. 

• PG 1211+143 

This host is too faint to be classified in this work. Taylor et al. (1996) have concluded 
from their ground-based if-band imaging that the host is dominated by a disk. 

• PG 1229+204 

This host was identified as a spiral from ground-based visible observations (Hutchings 
& Neff 1992) and HST/WFPC2 observations (Hamilton et al. 2002). Our observations 
show a prominent bar (radius is 5 kpc), with faint spiral arms starting from their ends. 
The result of the 2-D fit is a bulge+disk model about one magnitude brighter than 
the central nucleus. We classify this host as a "disk present" because of the structures 
visible in our image. 

• PG 1302-102 

The quasar is a strong flat-spectrum radio source. Two compact, but resolved, bright 
companions are within 10 kpc of the QSO. Although decoupling the closer companion 
and the host in the 2-D fitting analysis is difficult, visual inspection of the image 
suggests that the host of PG 1302-102 may be a somewhat asymmetrical elliptical 
galaxy. The 2-D fit favors a r 1 / 4 model about 1.5 magnitudes fainter than the central 
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nucleus, which we adopt for the classification of this host. Bahcall et al. (1997) and 
Hamilton et al. (2002) (working from the same //5T/WFPC2 image) both suggest an 
elliptical host. 

• PG 1307+085 

Our image shows a smooth symmetrical host with no apparent structure, and the 2- 
D fit shows that it is an elliptical galaxy with no detectable elongation. Analysis of 
HST/WFPC2 images gives a similar result (Bahcall et al. 1997; Hamilton et al. 2002). 
The four small ojects at the corners of a perfect 10" by 10" square centered on the 
QSO are sky subtraction artefacts. 

• PG 1309+355 

This host is smooth, symmetric and moderately elongated in the North-South direction. 
Although it has been described as an early-type spiral from HST/WFPC2 observations 
(Bahcall et al. 1997; Hamilton et al. 2002), our 2-D fit favors an elliptical host. HST 
observations also give a closely r 1,/4 -law profile, but pick up weak but definite spiral 
features, which are not seen in our image probably because of the longer wavelength. 
Classification for this host is therefore somewhat uncertain, although we adopt the 
"Elliptical" classification for consistency with the criteria adopted for the rest of our 
sample. 

• PG 1351+640 

Our image shows a large smooth host, slightly elongated in the north-south axis. 
Ground-based observations in the visible show a smooth host, best fitted as an el- 
liptical (Hutchings & Neff 1992). Our 2-D analysis does not allow classification of this 
galaxy. Deep optical imaging (Surace, Sanders, & Evans 2001) shows an faint exten- 
sion of the host to the west. This feature is marginally detected in our image, and 
appears as a subtle light excess about 5" west of the QSO. 

• PG 1411+442 

Highly disturbed host, with a bright elongated companion at a projected distance of 
3.6 kpc (2"1). Three long tidal arms or loops are visible in our image. This host had 
previously been classified as a spiral (Hutchings & Neff 1992). Deep optical images 
(Surace, Sanders, & Evans 2001) show that the tidal features extend out to at least 30 
kpc. Our 2-D fit does not allow classification of this host, and gives different answers 
in different bands. We therefore classify this object as highly disturbed. 

• PG 1426+015 

Strongly disturbed host, with a bright companion at a projected distance of 4.1 kpc 
(275). The host galaxy is asymmetrical: its south-west side is brighter and shows tidal 
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features. Our 2-D fit favors an elliptical host, but the apparent tidal features lead us 
to classify this object as strongly disturbed. 

• PG 1435-067 

Our image shows a smooth host galaxy. The northern half of the host is slightly 
brighter than the southern half at a radius of about 2 kpc. A companion is resolved 
(about 2 kpc diameter) at a projected distance of 10.8 kpc (4'.'8) to the west. Our 2-D 
fit suggests an elliptical, which we adopt for the host classification. 

• PG 1440+356 

Two symmetric spiral arms are visible in this host, without any obvious evidence of 
disturbance. The 2-D fit for this object is of poor quality, and its result (high surface 
brightness disk galaxy of small size) is probably driven by a mismatch in the PSF or 
the background sky level. Surace et al. (1998) classify the host as a spiral with a bar. 

• PG 1613+658 

This bright host galaxy shows a very complex structure, with prominent tidal arms and 
a bright resolved apparent companion projected 5.1 kpc (2'/2) west of the QSO. This 
structure was first reported by Hutchings & Neff (1992). Similar structure is shown 
by an HST/NICMOS image (McLeod & McLeod 2001). Deep optical images (Surace, 
Sanders, & Evans 2001) show tidal arms extending out to 35". 

• PG 1617+175 

The PSF subtraction is poor because of an APD failure in the AO system. Our image 
shows a smooth host and the apparent asymmetry in the image is most likely due 
to PSF subtraction residuals. An elliptical host is favored by the 2-D analysis, and 
adopted for host classification. Due to the poor quality of the PSF subtraction, this 
classification might not be reliable. 

• PG 1626+554 

Smooth host. The 2-D analysis suggests an elliptical host, but the reference PSF 
for this object is not very good, and the host classification might not be reliable. The 
rings visible in the residual wide-field image of Fig. ??g might be due to sky subtraction 
residuals. 

• PG 1700+518 

Highly disturbed host, with an interacting companion 7 kpc (1'.'6) north of the QSO. 
A tidal arc, especially bright in the north and west, is associated with the companion. 
Several bright knots are visible within this tidal feature. Fainter extended features 
to the south-west and south-east of the QSO are visible in our image. Despite the 
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high quality of our data, no reliable classification is obtained by our 2-D fit for this 
highly disturbed object. Stockton, Canalizo, & Close (1998) obtained imaging with 
an earlier version of the Hokupa'a AO system at the Canada-France-Hawaii telescope, 
detecting the companion and its arc. They also found spectroscopic evidence for a 
young (85 Myr) stellar population in the companion. Hines et al. (1999), from their 
iJST/NICMOS imaging of this object, suggest that the companion is a ring galaxy, 
produced by a small-impact-parameter collision. 

• PG 2130+099 

This spiral host shows two symmetric, high-contrast spiral arms extending out to at 
least 6'.'5 from the QSO. A non-stellar apparent companion is visible 9" to the east- 
southeast of the QSO. These features are also seen in ground-based and HST/WFPC2 
optical imaging (Hutchings & Neff 1992; Surace et al. 1998), which also show evidence 
for a central bulge. Strangely, Taylor et al. (1996) classify this object as an elliptical, 
although with marginal confidence. Despite the poor PSF subtraction in our image of 
this object, the 2-D fit correctly identifies this host contains a prominent disk. There 
is a good match between the fit parameters obtained in the H and K' images, both 
indicating that the host is almost 2 magnitudes fainter (<5m H = 1.56, 5m K > = 1.96) 
than the central nucleus. 

• 4C 31.63 

As the name indicates, this quasar is a strong radio source. The host is smooth and 
elongated, with 3 small apparent companions (projected distances are 11 kpc (2'.'5), 
12 kpc (2"8) and 20 kpc (4"7)) within 5" of the QSO. The companion 12 kpc to the 
South is of stellar appearance. Our image suggests that, in the central part of the 
host (at about 2" from the QSO), the southeast side of the galaxy is brighter than the 
northwest side. However, the effect is subtle and may be due to a PSF artifact; we 
therefore have chosen to classify this host as undisturbed. The 2-D analysis suggests 
that the host is an elliptical. The classification by Taylor et al. (1996) is ambiguous, 
whereas Hamilton et al. (2002) classify the host as an elliptical from an HST/WFPC2 
image. Spectroscopy of the host (Miller & Sheinis 2002) also support classification as 
an elliptical galaxy. 

• PG 2214+139 

This host has a very peculiar radial profile, although the 2-D analysis suggests that it 
is an elliptical. Sharp concentric rings surround the the central QSO at radii ranging 
from 2.5 kpc to 5 kpc, and are especially obvious in the north-west part of the host. 
These shells are probably tidal debris from a recent merger. 
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• PG 2251+113 

The quasar is a 4C and Parkes radio source. The smooth host is moderately elongated 
in the east-west direction, and it has a non-stellar apparent companion 4'.'3 to the 
North. The 2-D analysis shows a slight preference for a disk, but the classification 
is to be considered somewhat unreliable due to poor PSF subtraction. It would be 
quite unusual for a radio source of this power to have a host dominated by a disk, and 
Hutchings & Neff (1992) report an r 1 / 4 profile for it. 

• PG 2349-014 

This highly disturbed host shows 3 tidal arms and a close companion at 5 kpc (1'.'7) 
projected distance. The small companion is composed of a bright compact nucleus and 
two fainter knots about 0?!3 on either side of it, to the East and West. Spectroscopy 
(Miller & Sheinis 2002) has confirmed that the companion and QSO redshifts match. 
This host was identified as a elliptical undergoing strong interaction from //5T/WFPC2 
observations (Bahcall et al. 1997), which is in agreement with our 2-D fit and our 
image. Due to the features visible in our image, we chose to classify this host as 
strongly disturbed. The quasar itself is a strong Parkes radio source. 



4-4-3. Fit parameters for adopted host type 

Table 6 lists the value of the x 2 -model fit parameters (for the "adopted host" type listed 
in col. 8 of table 5) according for the following host decomposition : 

I(r, 6) = I e x io- 3 - 3307x ( - e ) + J x e ^ (2) 

where 

J(r, 9, e,PA)=rx^ 1.0 + sin*(9 - PA) x ^ j ^ - 1^ . (3) 

For all objects, we verified that the minimization algorithm gives results independent 
of the starting point, and that Xc xp+r i/4 is smaller than xj?i/4 an d Xcxp- For "E" and "D p " 
types in col. 8 of table 5 only the first or second term, respectively, on the right-hand side of 
equation 2 is relevant, while for type "B+D" all 9 parameters are used in the fit. When the 
host classification is uncertain ("?" in col. 8 of table 5), the 9 parameters of the bulge+disk 
fit decompositions are listed in Table 6. 



5. Discussion 



The results presented in the previous section illustrate both the overall success of our 
survey in being able to detect the great majority of the host galaxies of nearby PGQSOs, as 
well as the difficulties encountered when attempting to fit model radial surface brightness 
profiles to individual hosts. The two year timeframe of the survey provided sufficient oppor- 
tunity to revisit objects where the initial data had been compromised by poor seeing and/or 
problems with the stability of the PSF. In the end we were able to reliably detect the host 
galaxies in all but four of our targets (28/32), with sufficient signal-to-noise to allow fitting 
of the radial surface brightness profiles out to radii typically >10 kpc - a tribute in large part 
to the light gathering power of 8-m telescopes. However, despite extreme attention to detail, 
the stability of the PSF was still a limiting factor in preventing accurate subtraction of the 
central point source, thus compromising our ability to determine the host light profile a radii 
< 1-3 kpc, thus making the fit result rather sensitive to the choice of inner radius limit. Our 
model fits also proved to be somewhat sensitive to the exact choice of outer radius limit, 
thus in the end our refusal to use x 2_ fits alone, and our attempt to give at least equal weight 
to visual inspection of the images, in the final determination of "host type". Despite these 
caveats, there are still several important statistical results and trends that clearly emerge 
from the data. 



5.1. Classification statistics — host types and absolute magnitudes 

Twenty-eight of our 32 objects had detected hosts for which we attempted make a 
classification, and we discuss the distribution of "host type" and host absolute magnitude for 
these 28 objects below. The remaining 4 objects were not classified due either to having hosts 
which were too faint or having data which suffered from a bad PSF (see footnote to Table 7). 
Table 7 summarizes the number of objects and their mean H-band absolute magnitudes for 
each "host type" - "E", "D p ", "B+D" - as well as for those objects where a "host type" could 
not be reliably determined (i.e. "Unknown"), and also by "degree of disturbance" - "Non- 
disturbed", "Disturbed", "Strongly disturbed". These classifications follow the category 
listings as used in column "Dist." and column "Adopted type", respectively of Table 5. 

The percentage distributions of "host type" and the relationship between "degree of 
disturbance" and the "host type" can be seen from the graph in Fig. 6. It is immediately 
apparent that the hosts of the sample of PGQSOs in our survey are not characterized by 
a single dominant type. For the subsample of 19 objects with classified host type ("E", 
"Dp", "B+D" in table 5), there appears to be an almost equal split between objects with 
no detected disk component ("E") and objects with detected disks ("D p " and "B+D"). It 
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is also apparent that a significant fraction, (9/28 = 30%), of the hosts of PGQSOs show 
clear signs of disturbance and/or are strongly disturbed. If we then make the reasonable 
assumption that the "disturbed" , and even more likely the "strongly disturbed" , hosts also 
have a disk component, then the percentage of objects with a detectable disk component 
could be as large as 57% (16/28) to 68% (19/28). 

Table 7 also lists the median and average Mh values, and the la uncertainty in the 
average for both the host and nucleus, for the total sample (28 objects) and for each "host 
type" subsample. These values have been computed directly from the adopted fits for each 
object as listed previously in Table 6. The median and mean host magnitudes typically differ 
by less than 0.3 mag, with no obvious systematic offset. In what follows we will use the mean 
values to characterize the properties of the total sample and the "host type" subsamples. 

For our entire sample (28), the mean host absolute if-band magnitude is Mh = —24.82, 
which is equivalent to ~ 2.1 L H [assuming M H = 24.00 3 ]. Our value of ~ 2.1 L H for the mean 
host luminosity is almost identical to the values of ~ 2 and 2.3 found respectively by 
McLeod & Rieke (1994b), and Surace, Sanders, & Evans (2001) from their earlier ground- 
based observations of samples of QSOs selected from the same BQS sample of PGQSOs, 
and using the same z < 0.3 redshift cutoff as our current sample. The latter study also 
finds a range of 0.5-7.5 L H for the host galaxies, very similar to our measured range of 
-0.63-10.0 L H . 

Mean if-band host magnitudes and luminosities for samples of PGQSOs at z < 0.3, 
and for quasars in general, are not that abundant in the literature, but the few reported 
measurements all seem to fall generally in the range ~1-3L H , with the lower value being 
found for low-luminosity QSOs and the higher value being found for more luminous QSOs 
(e.g., McLeod & Rieke 1994a,b; McLeod & McLeod 2001). It is interesting to note that 
this range of host Mh is nearly identical to the range in the average Mh found for our 
different "host type" sub-samples. Table 7 shows that the lowest mean value — Mh = —23.78 
(~ 0.8L H ) — is found for hosts of type D p , i.e., hosts which have a dominant disk component, 
while the mean value for the ellipticals (type E) in our sample is M H = —25.24 (~ 3.1L H ). 
The similarity in range and mean //-band luminosities for different host type with previous 
mean values versus QSO luminosity could be explained, for example, if there were a clear 
correlation of host type with QSO luminosity, and indeed we do see such a correlation, as 
discussed below. 



3 Dcrivcd by adopting the mean observed value of H— K = 0.22 for the hosts and an observed value of M£ = 
—24.22 from Mobasher, Sharpies & Ellis (1993), corrected for our assumed value of H = 75kms _1 Mpc~ 1 . 
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5.2. Host properties vs. QSO luminosity 

A more complete representation of the range of properties of our sample of PGQSOs 
is given by the data shown in both Fig. 7 and Fig. 8. The color coding and shading of the 
symbols used to represent the individual objects is the same in each of the figure panels, and 
provides an aid to easily see the mean and range of host and nuclear properties for objects 
with different host type and degree of disturbance. 

Figure 7 (top panel) clearly shows a distinct change in host type with increasing M H of 
the host. Hosts with dominant disk components (i.e. "D p ") are preferentially clustered at 
the low end (i.e., ~sub-L H to L^) of the observed range, whereas E-type hosts dominate at 
above-L H . Objects which have both a relatively large bulge component and a detectable disk 
(i.e "B+D") are distributed more like E-type hosts, although the relatively small number 
(3) of "B+D" hosts limits our ability to say much more, other than that the they cluster 
near ~ L^. In general, Fig. 7 (top panel), shows that disturbed and strongly disturbed hosts 
are found at all values of host Mh, although there seems to be a tendency for the strongly 
disturbed hosts to be found at the upper end of our observed range of M B . Finally, we 
comment on the apparent bimodal distribution of host Mh in figure 7, simply by noting that 
with just a relatively small shift in the computed host Mh for a few objects, this could just 
as easily be a single distribution with a high luminosity tail. 

Figure 7 (bottom panel) also shows a change in the distribution of host type with 
increasing Mb of the QSO, similar to what is seen with increasing host Mh (top panel) — E- 
type hosts dominate at M B < —23 while at M B > —23 hosts with prominent disks ("D p " and 
"B+D") appear to become an increasing dominant fraction of all objects with decreasing 
B-band luminosity of the QSO. 

Several earlier imaging studies of low-redshift QSOs (i.e., z < 0.3) have also commented 
on the fact that lower luminosity QSOs appear to have a larger fraction of hosts with a 
distinguishable disk component, while high luminosity QSOs seem to reside primarily in 
elliptical hosts. The collective series of papers by McLeod and coworkers (McLeod & Rieke 
1994a,b, 1995; McLeod & McLeod 2001) used both ground-based and HST near-infrared 
imaging data of a PG sample that is most similar to ours to conclude that high luminosity 
QSOs seem to predominantly live in elliptical hosts while the host types of low-luminosity 
QSOs are more likely to show evidence for disks. Independently, Dunlop and coworkers 
(Taylor et al. 1996; McLure et al. 1999; Dunlop et al. 2003), using a combination of ground 
and space-based, near-infrared and optical imaging data for a mostly different sample of 
QSOs, reached generally similar conclusions when they emphasized the almost exclusive 
presence of elliptical hosts for the higher luminosity QSO in their sample, while also pointing 
out the increasing number of disks found in the hosts of the lower-luminosity, radio-quiet QSO 
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in their sample. And recently, Hamilton et al. (2002) provided a comprehensive summary of 
their HST archival study of 71 QSO host galaxies with WFPC2 images in F606W or redder 
filters, which concluded that "elliptical hosts are typically twice as luminous as spiral hosts" , 
although they also noted that the luminosity distributions of elliptical and spiral hosts for 
the radio-quiet QSOs were "generally compatible". 

A more comprehensive presentation of our data is given in Fig. 8, which shows the 
relationship of the computed 2-D fit value of Mh (column 5 of Table 6) and the observed 
Mb (column 5 of Table 1) for each QSO, along with the symbol for host type, the measured 
redshift and an indicator for radio loudness ("RLQ"). In addition to showing the trends 
of host type versus Mh and Mb discussed above, Fig. 8 also shows how the ratio of these 
two quantities (which was summarized by host type in cols. 9-11 of Table 7) varies over the 
sample. 

One interesting question is whether there is any evidence for a correlation between QSO 
luminosity and the luminosity of the host galaxy. The results shown in Fig. 8 suggest not, 
at least for hosts with —23.5 > M H (host) > —25.5 (i.e. ~ 0.6 - 4 L^(host)); if anything, 
the mean ratio M B /M H (host) versus M H (host), in this range of M H (host), seems to decrease 
with increasing Mr (host), but this result is not significant because of small number statistics. 
However, we also note that that all of our objects with Mn(host) > —25.5 are at z < 0.15, 
and the stated 5-band limit (B ~ 16.16) of the BQS (Schmidt & Green 1983) implies that 
we should still be complete out to this redshift for QSOs with Mb < —21.3. For hosts with 
Mn(host) < —25.5, Fig. 8 does seem to imply that they produce more luminous QSOs. 
However, this result is clearly compromised by incompleteness due to the 5-band limit of 
the BQS; all of the sources with M H (host) < —25.5 lie at redshift z > 0.13, which excludes 
detection of QSOs with M B > -23. 

As was pointed out by the referee, what does seem to be significant from Fig. 8 is that 
even the most luminous galaxies with disks do not seem to harbor really luminous QSOs. It 
would seem that it may be necessary to have a merger violent enough to completely destroy 
the disks of the progenitor galaxies to produce the most luminous QSOs. 

5.3. Radio-Loud vs. Radio Quiet hosts 

The number of radio-loud quasars (RLQs) in our sample is relatively small (5/28 = 18%), 
but it is representative of the fraction found in previous studies of large samples of QSOs, 
which show that RLQs are only ~15% of the total population of optically-selected QSOs. 
Although we did not specifically distinguish RLQs and RQQs when computing mean host H- 
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band luminosities for different "host types" in Table 7, the RLQ hosts are explicitly identified 
in Fig. 6-8. It is immediately obvious that RLQs live in the most luminous QSOs (i.e., 
M B < —24) with the most luminous hosts (i.e., M H (host) < —26 corresponding to > 6Lj^). 
These facts, coupled with our previously discussed trends of host type with QSO 5-band 
and host //-band luminosity, are consistent with previous published results which claim that 
RLQs are typically 1-2 mags more luminous than RQQs in the optical, and that RLQs are 
exclusively found in elliptical hosts, which in turn are usually much more luminous — typically 
1-2 magnitudes, depending on whether the observations are in the optical or near-infrared, 
respectively — than the mean values found for the hosts of RQQs (e.g., Hutchings, Crampton, 
& Campbell 1984; Smith et al. 1986; McLure et al. 1999; Hamilton et al. 2002; Dunlop et al. 
2003). 

Although the origin of the observed dichotomy in the mean properties of RLQ and RQQ 
hosts is somewhat beyond the scope of our current study, it is nevertheless worth restating the 
obvious - that QSO surveys biased toward radio-selected objects will necessarily find more 
luminous, elliptical hosts when compared with pure optically selected samples of QSOs. 

5.4. Degree of disturbance of host galaxy 

One of the expected strengths of AO observations of QSOs using 8m telescopes was 
the ability to better characterize the "host type" , and in addition to detect faint extended 
structure in the host which may have escaped detection in previous imaging observations with 
smaller telescopes. Although the majority of our host classifications are in general agreement 
with previous results (see §4.4.2), our images were more likely to allow us to distinguish 
asymmetric features (e.g. putative tidal debris, central bars, etc.) from, for example, more 
symmetric spiral structure, and to better separate out compact, circumnuclear structures 
from the bright QSO nucleus. We note that since we rely on morphological features such 
as tidal arms to classify hosts as "not disturbed" , "disturbed" or "strongly disturbed" , this 
classification is relatively insensitive to usually very symmetrical PSF subtraction residuals. 
For example, the host type of PG 1626+554 and PG 2251+113 are unknown due to poor 
PSF subtraction, but the smooth host images obtained lead us to classify them as "not 
disturbed" . 

According to our analysis, 30% (9/30) of the classified QSO hosts were categorized as 
being "disturbed" or "strongly disturbed" in Table 5. The 6 hosts classified as "strongly 
disturbed" show asymmetric features reminiscent of tidal structure, as might be due, for 
example, to the advanced merger of two spiral disks. Putative "tidal arms" are clearly seen 
in PG0157+001, PG 1411+442, PG 1613+658, PG 1700+518 and PG 2349-014, while PG 
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1426+215 has a clearly asymmetric disk. The three additional objects listed as "disturbed" 
generally show fainter asymmetric features, (plausibly faint tidal debris) whose effect on the 
1-D radial profile is small enough to still allow us to distinguish between different host types. 
The 4 "disturbed" hosts include two "D p " and one "E". 

In addition to the 9 hosts classified as either "disturbed" or "strongly disturbed" , there 
are 3 hosts that were found to have a prominent well-defined bar, whose size in each case was 
~ 5 kpc radius. In two of these cases the bar actually dominated the integrated luminosity 
of the host. The hosts of two of these barred galaxies (PG0804+761, PG1001+054), were 
classified as having a dominant disk profile ("D p "), while the third (PG1229+204), was 
classified as "B+D" . It is unclear whether these bars have any relation to the processes that 
produced the putative tidal features in the "disturbed" and "strongly disturbed" hosts, but 
if they do, then the percentage of such hosts rises to 40% (12/30). 

Our finding of 30% for the fraction of "strongly disturbed" + "disturbed" hosts (40% if 
"well-defined bars" are included) is similar to the range of values, typically 25%-50%, quoted 
for "morphologically peculiar", and/or "interacting" hosts in previous surveys of low-z QSOs 
(e.g., Hutchings, Crampton, & Campbell 1984; Smith et al. 1986; McLeod & Rieke 1994b; 
Marquez et al. 2001; Surace, Sanders, & Evans 2001). Although it is not easy to determine 
if the definitions used to define degree of disturbance in previous surveys are the same as 
used here, a comparison of those objects from our current survey which have been included 
in previous surveys suggests that we are in fact using similar criteria. 

The breakdown by "host type" of "disturbed" , "strongly disturbed" and "well defined 
bar" systems is given in Fig. 6. Only one such object (PG 2214+139) with obvious "concen- 
tric shells" is classified as an elliptical. The remainder of the "disturbed" and "well defined 
bar" systems have obvious disks, as presumably do the "strongly disturbed" objects (see 
above). Taken together, 11 of these 17 systems (65%) have features which are interpreted 
as being indicative of obviously disturbed, non-axisymmetric structure. Such structure has 
in the past been seen as evidence for interaction/mergers as a possible trigger for producing 
the QSO. 

The breakdown by M H (host) and M B (total) of "disturbed", "strongly disturbed" and 
"well defined bar" systems is given in Fig. 7. In sum, there appears to be no clear evidence 
for a correlation with either M H (host) or M B (total) when these three classes are considered 
together. However, if we consider only the "strongly disturbed" systems, there is a clear 
preference for these to be found in the most luminous (i.e., > 2L^) hosts. 
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5.4-1- Close compact companions or super- star- cluster knots 

In addition to asymmetric structures, we also found a small, bright knot in the inner 
disks of 9 hosts (see notes in Table 5). These knots were typically within a few kpc (projected 
distance) of the QSO. They are being discussed here since we believe that the evidence points 
to the knots being associated directly with the host galaxy, as opposed, for example, to being 
projected foreground or background objects at much larger distances from the QSO. 

The terms "companion" or "close companion" have been widely used in the literature 
to describe a range of distance scales around QSOs where galaxies at the same redshift as 
the QSO are found. For example comparative studies of the environment around QSOs 
typically find that, at low redshifts, QSOs live in regions of moderate galaxy overdensity, 
typical of compact or small groups (e.g., Yee & Green 1984; Heckman et al. 1984). Closer in, 
typically within a projected nuclear separation of < 30 — 50 kpc between the companion and 
QSO nuclei, are numerous reports of close companion galaxies, clearly distinguishable as a 
separate galaxy, but often with reported evidence of a tidal interaction with the QSO host 
(e.g., Stockton 1982; Heckman et al. 1984; Hutchings, Crampton, & Campbell 1984; Bahcall 
et al. 1997; Boyce, Disney, & Bleaken 1999; Marquez et al. 2001). At the closest distances 
(typically nuclear separations < 10 kpc), where the putative close companion galaxy strongly 
overlaps the QSO host, as is the case in most of the strongly disturbed systems, including 
our own, the close companion is then plausibly the second object involved in an advanced 
merger with the QSO host. This last case is similar to what is meant here by the term "close 
companion" , except in our case we also clearly distinguish a relatively bright "compact knot" 
which could either be the nucleus of the companion galaxy, or possibly a super star-cluster, 
perhaps triggered by the merger. 

We have referred to these objects (the ones within lOkpc projected separation) as "close 
compact companions" even though redshift is available for only one of the 9 "knots", since 
it is likely that these "knots" are indeed at the same redshift as the QSO rather than simply 
being unrelated foreground or background objects. Five of the 9 QSO hosts with "close 
compact companions" are the same 5 hosts classified as being "strongly disturbed" , and one 
of the remaining four is classified as being "disturbed" . Several of the companions are also 
observed to be strongly elongated (e.g., in PG 1119+120, PG 1411+441, PG 1426+215, PG 
2349—014). These facts strongly suggest that the knots are directly related to a merger 
event. 

However we remain somewhat agnostic whether these objects represent the second nu- 
cleus of a nearly merged galaxy pair, or, instead, a super star cluster perhaps triggered by 
the merger. Although, the "knots" appear to be too bright to be "normal" star-formation 
knots, and are therefore good candidates to be distinct galaxy nuclei, they are in fact sim- 
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ilar in brightness to the putative "super star cluster" knots that have often been found in 
the circumnuclear regions of "infrared-excess QSOs" (e.g., Surace, Sanders, & Evans 2001) 
and that appear to be nearly ubiquitous in the inner disks of "warm ultraluminous infrared 
galaxies" (e.g. Surace & Sanders 1999). All but one of the "knots" seen in our current QSO 
sample have measured Mr in the range —20.1 to —22.3, and in each case this translates into 
an absolute if -band luminosity in the range 3.5-5.5 mag fainter than the host galaxy, and 
the most luminous knot (M H = —23.18), is still 3.3 mag fainter than its host (PG1302— 102), 
and is still within the upper range of luminosities observed for the star forming knots in 
"warm ULIGs". 



5.5. Mid/Far-infrared properties of QSO hosts 

In addition to trying to better categorize the host types for our sample of QSOs, we 
were also interested in investigating their mid- and far-infrared emission properties, and in 
particular to see if a relationship exists between the infrared luminosity of the QSO and 
morphology of the host. 

Interest in the infrared properties of QSOs has been fueled by suggestions of an evo- 
lutionary connection between QSOs and ultraluminous infrared galaxies (ULIRGs), where 
both types of objects represent different phases in the end stage of the merger of two rela- 
tively equal mass gas-rich spirals (see review by Sanders & Mirabel 1996). Key to such a 
connection has been the suggestion that those QSOs with the strongest mid/far- "infrared 
excess" emission seem to have "disturbed" hosts (e.g., Sanders & Mirabel 1996; Boyce et 
al. 1996; Canalizo & Stockton 2000; Clements 2000; Haas et al. 2000; Surace, Sanders, & 
Evans 2001), whose properties are somewhat similar to the hosts of "warm" ULIRGs. For 
the QSOs, "infrared-excess" has been defined using the mean radio-to- Xray spectral en- 
ergy distributions (SEDs) of QSOs (Sanders et al. 1989; Elvis et al. 1994). The data have 
been interpreted as showing that the SED is dominated by two thermal emission "bumps" . 
In addition to the well-known "blue bump" (BB) of optical/UV emission at wavelengths 
~ 0.01 — 1.0 //m, thought to be associated with thermal emission from an accretion disk, 
there also exists an "infrared/submillimeter bump" at wavelengths ~ 1.0 — 500 /im, which is 
interpreted as being primarily due to dust re-radiation of emission from either the accretion 
disk or from embedded circumnuclear star formation. The "infrared excess" is then defined 
to be the ratio of luminosities in the two "bumps", Lir/Lbb- 

The recent availability of reprocessed ISO observations of QSOs (Haas et al. 2000; 
Polletta et al. 2000; Haas et al. 2003), plus the release of the IRAS Faint Source Catalog 
(Moshir et al. 1992), have now made it possible to more accurately estimate the strength as 
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well as the color temperature of the "infrared/submillimeter bump" for the BQS QSOs. In 
particular, we have been able for the first time to compile fluxes in four wavebands (12/im, 
25/iin, 60/zm and 100/im) for the majority (18/32) of our BQS sources, and in 2-3 of these 
wavebands for an additional 10 sources. Three sources have only one flux measurement, but 
with sensitive upper limits in the remaining 3 wavebands. Only one source has just upper 
limits in all four wavebands. These measurements are listed here in Table 8 where they have 
been used to compute an infrared luminosity, L IR = L(8 — 1000)//m, and the luminosity 
ratio, Li R /L BB , for each source following the definitions and prescriptions found in Table 1 
of Sanders & Mirabel (1996). 

Using the data in Table 8, we compare the infrared/optical ratio, L IR /L BB , of the QSO 
hosts versus host type in Fig. 9, and tabulate the mean and median values of L IR and median 
Lir/Lbb, for different host type and degree of disturbance in Table 9. Figure 9 suggests 
a fairly strong correlation between Li R /L B b and host type, with nearly all ellipticals ("E") 
having the lowest observed values of Lui/ Lbb- Table 9 shows that Es have median values of 
Ltr/Lbb slightly less than half that of hosts with detectable disks ("D p " + "B+D"). The 
previously observed correlation between non-disturbed hosts and E host type then implies (as 
shown in Table 9) that the median L^/Lbb ratio for "disturbed" and "strongly disturbed" 
hosts is also approximately twice the median ratio for "non-disturbed" hosts. 

A strong correlation between IR luminosity and host morphology had previously been 
suggested by Clements (2000), although a close comparison between his "quiet/disturbed" 
classification and the images presented here shows that the images he used lacked the angular 
resolution and/or depth to reliably classify the host galaxies: for example, PG1411+442 and 
PG1700+518, two of the most disturbed objects in our sample, were classified as "quiet", 
while PG0923+201 and PG0844+349 were classified as disturbed. 

Another method of characterizing the mid/ far- infrared emission of QSOs has been to use 
the mid/far- infrared color temperature (e.g., Lipari 1994; Canalizo & Stockton 2001). The 
color-color plot shown in Fig. 10 adopts the same definitions used in Canalizo & Stockton 
(2001), who first pointed out that nearly all (8/9) of the "transition QSO" they studied — i.e., 
objects with an "intermediate position" in the far-infrared color-color diagram, in-between 
the empirical location of optically-selected QSOs and far-infrared-selected ULIRGs — were 
undergoing strong tidal interactions indicative of major mergers. Our new results are similar, 
in that all 8 of our "disturbed" and "strongly disturbed" QSOs for which we could reliably 
compute far-infrared colors, fall in the "transition region" between the two limits marked by 
the lines denoting "power-law" and "black-body" emission. Additionally, we find only objects 
we classified as "non-disturbed" above the "power-law" continuum line. We interpret this 
as additional evidence for a positive correlation between fading evidence of tidal structure 
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(i.e., increasing age of the merger), and decreasing Lir/LbB) since the latter is correlated 
with increasing far-infrared color temperature. 

6. Conclusions 

We have used AO imaging on the 8.1m Gemini-North and 8.2m Subaru telescopes to 
acquire near-infrared (mostly if-band) images of an unbiased sample of 32 nearby (z < 
0.3) optically selected QSOs from the PGBQS in order to investigate the luminosity and 
morphology of the host galaxies, and to compare the properties of the host with properties 
of the QSO. The 5-band absolute magnitudes of our QSO targets span a factor of ~ 23, from 
~ Mb — —22.12 (the minimum threshold for "bona-fide QSOs" according to the definition 
given by Schmidt & Green 1983) to Mb = —25.53. Three sources have absolute luminosities 
slightly below the minimum threshold and are traditionally classified as "Seyfert 1 nuclei" . 
Given that this threshold is somewhat arbitrary, we have included these three objects in 
the current study. Following the definitions of Kellermann et al. (1994), twenty-seven of our 
sources were radio-quiet and 5 were radio-loud. 

Our imaging campaign in general was quite successful, given that we were able to detect 
and characterize the host type for 28 of our sources. Our basic results can be summarized 
as follows: 

1. The mean if-band magnitude (and corresponding L^) of the host galaxies for our 
sample of QSOs was M H (host) = -24.82 (~ 2.1 L^), with a range of -23.5 to -26.5 
(~ 0.63- 10.0 L* H ). 

2. For the sample of 28 classified hosts, the distribution of "best-fit" host types was 10 
(36%) ellipticals ("E"), 3 (11%) bulge+disk ("B+D"), 8 (29%) disk dominated ("D p "), 
and 7 (25%) of indeterminate type. 

3. There appears to be a strong correlation between host type and the iJ-band absolute 
magnitude of the host. Sub-LJj hosts all have a dominant ("D p ") or strong disk com- 
ponent ("B+D"); 1 — 2Ly hosts appear to be equally divided between ellipticals ("E") 
and disks ("D p " or "B+D"); > 2L^ hosts are mostly ellipticals ("E"). A similar, but 
somewhat weaker trend is found for the £>-band absolute magnitude of the QSO. 

4. Thirty percent (9/30) of our detected hosts were classified as either being "disturbed" 
(3) or "strongly disturbed" (6). An additional 3 objects were found to have strong 
well defined bars, bringing the total fraction of objects with obvious disturbed and/or 
non-axisymmetric structure up to 40% (12/30). In total, these objects are relatively 
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equally distributed over all values of M H (host), and M B (total), although in detail, 
the "strongly disturbed" objects are much more likely to be found among the most 
luminous hosts (i.e. > 2LJj), with "disturbed" and "well defined bars" more likely to 
be associated with less luminous hosts (i.e. < 2L^). 

5. For the subsample of 18 QSOs with redshifts z < 0.13, which is the redshift com- 
pleteness limit corresponding to the m Blim ~ 16.2 limit of the PGBQS (Schmidt & 
Green 1983), we find no obvious correlation between the absolute magnitude of the 
QSO, MB(total), and the absolute magnitude of the host galaxy, Mn(host) . How- 
ever at z > 0.15, all 8 sources in our sample have the largest host luminosities, 
M H (host) < -25.5 ( > 4I£), and the largest QSO luminosities, M B (total) < -23.8, 
suggesting that the most luminous QSOs are found among the most luminous hosts. 

6. The mean properties of the hosts of radio-loud quasars (RLQs) are clearly different 
from the mean properties of the hosts of radio-quiet QSOs (RQQs). All 5 of the RLQs 
are among the 10 sources in our sample at z > 0.17. All 5 RLQs are among the 10 
most luminous QSOs (i.e. Mb < —24), and for the 4 RLQs where we were able to 
reliable measure the host, all 4 are among the 5 most luminous hosts (i.e. > 6L^). 
Three of the 5 RLQ hosts are classified as ellipticals, one host suffers from a bad PSF 
but otherwise shows no visible sign of disk structure, and one is strongly disturbed and 
cannot be classified. These results are in keeping with the paradigm which suggest that 
RLQs are preferentially found in giant elliptical hosts with the most luminous QSOs. 

7. We find a strong correlation between the infrared excess, Lir/Lbb, of the QSO and its 
host type and host degree of disturbance. This ratio is twice as large, on average, in 
hosts with strong disk components as in elliptical hosts, and likewise is twice as large, 
on average, in "disturbed" + "strongly disturbed" hosts as in "non-disturbed" hosts. 
We also find that "disturbed" and "strongly disturbed" hosts have mid/far-infrared 
colors that place them in a region of the far-infrared color-color plane in between the 
power-law that is characteristic of the majority of QSOs, and the black-body that is 
characteristic of the majority of ULIRGs. 

Our results concerning the mean and range of host luminosities, the general breakdown 
of host types, the fraction of objects with disturbed features, and the differences in the mean 
properties of the hosts of RLQs versus RQQs, are in general agreement with previous studies 
of the host properties of low-z QSOs. Our mean and range of if-band host luminosities are 
almost identical to the values reported earlier by McLeod & Rieke (1994b) and Surace & 
Sanders (2000), respectively, in their studies of a similar sample of PGQSOs selected from 
the BQS (Schmidt & Green 1983). One advantage of our AO observations with 8m-class 
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telescopes is that we have perhaps been able to give more emphasis to trying to accurately 
characterize the host properties of individual sources, and to search for systematic trends in 
host type versus other properties of the QSO. 

Perhaps the most interesting new results from our survey are those concerning the 
detailed distribution of disturbed hosts versus the luminosity of the host and the luminosity 
of the QSO, and the apparent correlation of infrared excess with host type. Our results 
suggest that the percentage of disturbed hosts is significant (15-30%) across the range of 
observed host luminosities (~ 0.6 — 10 L^), even though the main host type clearly changes 
from hosts with prominent disks to hosts which appear to be luminous ellipticals. Whether 
disturbed hosts represent a particular phase is the lifetime of virtually all QSOs or whether 
they represent an inherently separate class of QSOs is somewhat unclear. The answer to 
this question depends on the answers to two subsidiary questions: (1) What is the time scale 
for QSO activity relative to that for the persistence of recognizable signs of disturbance due 
to a strong interaction or merger? and (2) What is the true space density of QSOs with 
bolometric luminosities above some specific value and the fraction of these that show some 
level of disturbance indicating a recent merger? Our study has some bearing on this second 
question. However, it now seems likely that UV-selected optical QSO samples, such as the 
one studied here, may represent considerably less than half of the total QSO population at a 
comparable bolometric luminosity. In addition, recent work by Jester et al. (2005) comparing 
QSOs found in the Sloan Digital Sky Survey with the PG sample shows that the PG sample is 
seriously incomplete even in terms of its own specified selection criteria. This incompleteness 
suggests two imperatives for future work. First, it seems essential to develop a sample that is 
substantially complete to a specific bolometric QSO luminosity over some large volume. This 
effort will necessarily require complementary techniques covering wavelength regimes at least 
from X-rays to the far IR. Second, it will be necessary to carry out imaging studies similar 
to those we have done for the PGBQS on representative subsamples of QSOs identified by 
these other techniques, for it is quite likely that, for example, the fraction showing signs of 
mergers may vary for samples selected by different criteria. A recent step in this direction 
has been taken by Hutchings et al. (2003) for QSOs with z < 0.3 found in the 2MASS survey, 
for which they find an extraordinarily high fraction (> 70%) show disturbances indicating 
tidal interactions. 

We thank the referee for a careful reading of the paper and useful comments to help 
improve both its content and its presentation. We also thank Gabriela Canalizo for a helpful 
conversation. The Hokupa'a AO observations were supported by members of the University 
of Hawaii AO group: P. Baudoz, O. Guyon and D. Potter. Support for Hokupa'a came from 
the National Science Foundation under grant no. AST 96-18852. Our results were based in 
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of Japan. 

REFERENCES 

Andreani, P., Franceschini, A., & Granato, G. 1999, MNRAS, 306, 161 
Bahcall, J.N., Kirhakos, S., Saxe, D. H., & Schneider, D. P. 1997, ApJ, 479, 642 
Balzano, V. A., & Weedman, D. W., ApJ, 243, 756 

Barthel, P. B., & Sanders, D. B. eds., 2006, QSO Host Galaxies: Evolution and Environment, 
New A Rev., in press 

Bothun, G. D., Romanishin, W., Strom, S. E., & Strom, K. M. 1984, AJ, 89, 1300 

Boyce, P. J., Disney, M. J., & Bleaken, D. G. 1999, MNRAS, 302, L39 

Boyce, P. J. et al. 1996, ApJ, 473, 760 

Canalizo, G., & Stockton, A. 2000, A J, 120, 1750 

Canalizo, G., & Stockton, A. 2001, ApJ, 555, 719 

Clements, D.L. 2000, MNRAS, 311, 833 

Clements, D. L., & Perez-Fournon, eds. 2001, Quasar Hosts, (Berlin: Springer), 336pp 

Cutri, R. M., Nelson, B. O., Francis, P. J., & Smith, P. S. 2002, in AGN Surveys, eds. R. 
Green, E. Khachikian, D. Sanders, (San Francisco: ASP Conf. Ser.), 284, 127 

Disney, M. J., et al. 1995, Nature, 376, 150 

Dunlop, J. S., Taylor, G. L., Hughes, D. H., & Robson, E. I. 1993, MNRAS, 264, 455 

Dunlop, J. S., McLure, R. J., Kukula, M. J., Baum, S. A., O'Dea, C. P., & Hughes, D. H. 
2003, MNRAS, 340, 1095 



Elvis, M. et al. 1994, ApJS, 95, 1 

Giveon, U., Maoz, D., Kaspi, S., Netzer, H., & Smith, P. S. 1999, MNRAS, 306, 637 

Graves, J. E., Northcott, M. J., Roddier, F. J., Roddier, C. A., Close, L. M. 1998, Proc. SPIE, 
3353, 34 

Haas, M., Klaas, U., Coulson, L, Thommes, E., & Xu, C. 2000, A&A, 356, 83 
Haas, M., et al. 2003, A&A, 402, 87 

Hamilton, T. S., Casertano, S., & Turnshek, D. A. 2002, ApJ, 576, 61 

Heckman, T. M., Bothun, G. D., Balick, B., & Smith, E. P. 1984, AJ, 89, 958 

Hines, D. C, Low, F. J., Thompson, R. I., Weymann, R. J., & Storrie-Lombardi, L. J. 1999, 
ApJ, 512, 140 

Hutchings, J. B., Crampton, D., & Campbell, B. 1984, ApJ, 280, 41 
Hutchings, J. B., Janson, T., & Neff, S. G. 1989, ApJ, 342, 660 
Hutchings, J. B., & Crampton, D. 1990, AJ, 99, 37 
Hutchings, J. B., & Neff, S. G. 1991, AJ, 101, 434 
Hutchings, J. B., & Neff, S. G. 1992, AJ, 104, 1 
Hutchings, J. B., & Neff, S. G. 1996, AJ, 113, 550 

Hutchings, J. B., Maddox, S., Cutri, R. M., & Nelson, B. O. 2003, AJ, 126, 63 
Hyland, A. R., & Allen, D. A. 1982, MNRAS, 199, 943 
Jahnke, K., & Wisotzki, L. 2003, MNRAS, 346, 304 
Jester, S., et al. 2005, AJ, 130, 873 

Kellermann, K. I., Sramek, R. A., Schmidt, M., Green, R. F., & Shaffer, D. B. 1994, AJ, 
108, 1163 

Lipari, S. 1994, ApJ, 436, 102 

Low, F. J., Cutri, R. M., Huchra, J. P., & Kleinmann, S. G. 1988, ApJ, 327, L41 
MacKenty, J. W., & Stockton, A. 1984, ApJ, 283, 64 



-38- 

Marquez, L, Masegosa, J., del Olmo, A., Lara, L., Garcia, E., & Molina, J. eds., 2001, QSO 
Hosts and Their Environments, (New York: Kluwer), 379pp 

Marquez, L, Petitjean, P., Theodore, B., Bremer, M, Monnet, G., & Beuzit, J.-L., A&A, 
371, 97 

McAlary, C. W., McLaren, R. A., McGonegal, & R. J., Maza, J. 1983, ApJS, 52, 341 
McLeod K. K., & McLeod B. A. 2001, ApJ, 546, 782 
McLeod, K. K., & Rieke, G. H. 1994a, ApJ, 420, 58 
McLeod, K. K., & Rieke, G. H. 1994b, ApJ, 431, 137 
McLeod, K. K., & Rieke, G. H. 1995, ApJ, 454, L77 

McLeod, K. K., Rieke, G. H., & Storrie-Lombardi, L. J. 1999, ApJ, 511, L67 

McLure, R. J., Kukula, M. J., Dunlop, J. S., Baum, S. A., O'Dea, C. P., & Hughes, D. H. 
1999, MNRAS, 308, 377 

Miller, J. S., & Sheinis, A. I. 2002 , ApJ, 588, L9 

Mobasher, B., Sharpies, R. M., & Ellis, R. S. 1993, MNRAS, 263, 560 

Moshir, M., et al. 1992, Explanatory Supplement to the IRAS Faint Source Catalog, Version 
2 (JPL D-10015) (Pasadena: JPL) 

Neugebauer, G., Oke, J. B., Becklin, E. E., & Matthews, K. 1979, ApJ, 230, 79 

Neugebauer, G., Green, R. F., Matthews, K., Schmidt, M., Soifer, B. T., & Bennett, J. 1987, 
ApJS, 63, 615 

Neugebauer, G., Soifer, B. T., Matthews, K., & Elias, J.H. 1989, ApJ, 97, 957 

Percival, W. J., Miller, L., McLure, R. J., & Dunlop, J. S. 2001, MNRAS, 322, 843 

Polletta, M., Courvoisier, T. J.-L., Hooper, E. J., & Wilkes, B. J. 2000, A&A, 362, 75 

Rudy, R. J., Levan, P. D., & Rodriguez-Espinosa, J. M. 1982, A J, 87, 598 

Sanders, D. B., Soifer, B. T., Elias, J. H., Madore, B. F., Matthews, K., Neugebauer, G., & 
Scoville, N. Z. 1988, ApJ, 325, 74 

Sanders, D. B., Phinney, E. S., Neugebauer, G., Soifer, B. T., & Matthews, K. 1989, ApJ, 
347, 29 



-39- 



Sanders, D. B., & Mirabel, I. F. 1996, ARA&A, 34, 749 
Schmidt, M., & Green, R. F. 1983, ApJ, 269, 352 

Scoville, N. Z., Frayer, D. T., Schinnerer, E., & Christopher, M. 2003, ApJ, 585, L105 

Sitko, M. L., Stein, W. A., Zhang, Y.-X., & Wisniewski, W. Z. 1982, ApJ, 259, 486 

Smith, E. P., Heckman, T. M., Bothun, G. D., Romanishin, W., & Balick, B. 1986, ApJ, 
306, 64 

Staguhn, J. G., Schinnerer, E., Eckart, A., k Scharwachter, J. 2004, ApJ, 609, 85 
Stockton, A. 1982, ApJ, 257, 33 

Stockton, A., Canalizo, G., & Close, L. 1998, ApJ, 500, L121 
Stockton, A., & MacKenty, J. W. 1983, Nature, 305, 678 
Stockton, A., & MacKenty, J. W. 1987, ApJ, 316, 584 

Surace, J. A., Sanders, D. B., Vacca, W. D., Veilleux, S., & Mazzarella, J. M. 1998, ApJ, 
492, 116 

Surace, J. A., & Sanders, D. B. 1999, ApJ, 512, 162 

Surace, J. A., & Sanders, D. B. 2000, AJ, 120, 604 

Surace, J. A., Sanders, D. B., & Evans, A. S. 2001, AJ, 122, 2791 

Takami, H., et al. 2004, PAS J, 56, 225 

Taylor, G. L., Dunlop, J. S., Hughes, D. H., & Robson, E. I. 1996, MNRAS, 283, 930 
Tokunaga, A. T., Simons, D. A., & Vacca, W. D. 2002, PASP, 114, 180 
Veilleux, S., et al. 2005, ApJ, submitted 
Veron-Cetty, M.-P, & Woltjer, L. 1990, A&A, 236, 69 
Yee, H. K. C, & Green, R. F. 1984, ApJ, 280, 79 

Zheng, Z., Wu, H., Mao, S., Xia, X.-Y., Deng, Z.-G., & Zou, Z.-L. 1999, A&A, 349, 735 



This preprint was prepared with the AAS IATj^X macros v5.2. 



-40- 




1e-06 



100 



Radius in pixels (1 pixel = 0.02") 



Fig. 1. — PSF radial profile variations. The 2 radial profiles plotted in this figure (PSF 1 and 
PSF 2) are derived from 30 s exposures of the same star taken 5 min apart and demonstrate 
the changes in the PSF radial profile due to seeing variations. In the small inset, the log- log 
radial profile of 2 PSFs of different brightness (PSF 4 is brighter than PSF 3) but under the 
same seeing conditions illustrate the effect of the PSF brightness on the radial profile. 
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Fig. 2. — Effect of the Cassegrain rotator angle variations on the PSF stability. The left 
image is the difference between two consecutive 10 s images of a star without any Cassegrain 
rotator rotation. The right image is the difference between two consecutive 10 s images 
between which the Cassegrain rotator was rotated by 15°. The intensity scale is identical in 
both images. 
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Fig. 3. — Result of the coaddition of the 33 frames (30 s exposure per frame) of PG1700+518 
(left). On the right, the total exposure for each pixel of the image is represented and ranges 
from 990 s (central part) to s (corners). 



-43- 




Fig. 4. — Raw (left), PSF-subtracted (center) and PSF+host model subtracted (right) images 
of the QSOs. The images are 10" x 10" (left and center) and 20" x 20" (right). The black 
bar in the left images is 5 kpc long. North is up, East is left on all images. 




Fig. 4. — Continued. 
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Fig. 4. — Continued. 
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Continued. 
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Continued. 




Fig. 4. — Continued. 




Fig. 4. — Continued. 
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Fig. 5. — 1-D radial profiles of the 2-D raw images and reference PSF (top), and residual 
errors of the r 1//4 and exponential 2-D fits. 
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Fig. 5. — Continued. 
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Fig. 6. — Distribution of "host type", and "degree of disturbance" within each "host type" 
for 28 out of the 32 objects (28 = 32 - 2 non detections - 2 "bad PSF") listed in Table 5. 
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Fig. 7. — (top) Distribution of host types as a function of the absolute if -band magnitude 
of the host. Only those 27 objects with hosts characterized (32 images — 2 "too faint to 
measure host" — 2 "bad PSF" — PG0050+124) are included in the histogram. If no H band 
images were acquired, H magnitudes for the QSO and host galaxy were derived from J or 
K', assuming J — H = 0.94 and H — K = 0.97 for the QSO nucleus (average values derived 
from the photometry obtained in our sample). For the host galaxy, J — H = 0.72 — 0.4 x z 
and H — K = 0.22 + 2.0 x z were assumed, where z is the redshift. (bottom) Distribution 
of host types as a function of the absolute total i?-band magnitude for all 32 QSOs. As 
described in the text, host types (disk present, bulge+disk or elliptical) are not assigned to 
strongly disturbed hosts. 
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Fig. 8. — Absolute total -B-band magnitude of the QSO versus the absolute if -band magni- 
tude of the QSO host galaxy. Only those 27 objects with hosts detected and characterized 
(32 images - 2 "too faint to measure host" - 2 "bad PSF" - 1 "no reliable 2-D fit") are 
shown (same subsample as top panel in Fig. 7). If no H band images were acquired, H 
magnitudes for the QSO and host galaxy were derived from J or K', assuming J—H = 0.94 
and H—K = 0.97 for the QSO nucleus (average values derived from the photometry obtained 
in our sample). For the host galaxy, J—H = 0.72 — 0.4 x z and H—K = 0.22 + 2.0 x z were 
assumed, where z is the redshift. Radio-loud objects are labeled in bold red. 
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Fig. 9. — Distribution of host types as a function of "infrared" (Lir) to "optical/UV" (Lbb) 
luminosity ratio of the QSO, where Lir and Lbb are meant to represent the total far-infrared 
luminosity and the luminosity of the "blue-bump" in the spectral energy distributions of 
QSOs (see text for the exact frequency limits used to estimate Lir and Lbb)- As described 
in the text, host types (disk present, bulge+disk or elliptical) are not assigned to strongly 
disturbed hosts. All 32 objects are shown in this figure, and PG2251+113 (Li R /L B b < 0.25) 
is shown here in the lowest L ir /L B b bin. 
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Fig. 10. — IR color-color diagram (adapted from Lipari (1994) and Canalizo & Stock- 
ton (2001)). Spectral indexes a(100,60) = — log(/ino M m//60/xm)/log(100/im/60/im) and 
a(60,25) = — log(/60/im//25/im)/log(60/xm/25/xm) are derived from Table 8. The light blue 
and light red areas show the approximate empirical locations of QSOs (close to a power law 
SED) and ULIRGs (dominated by thermal black body emission peaking in the mid-IR). Due 
to incomplete IR photometry, only 23 out of 32 objects are shown in this figure. Objects not 
detected (upper limit) in only one of the 3 IR wavelengths are shown, with a black arrow 
indicating possible positions on this figure. Lines corresponding to power-law emission and 
single temperature black body emission are also shown. Radio-loud objects are labeled in 
bold red. 
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Table 1. Source List 
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.4 


± 





,l d 


-25. 


.35 


Q 


PG2130+099 


21 


32 


27 


.92 


+ 10 


08 


18, 


.7 


0. 


,063 


14 


.7 


± 


0. 


l d 


-22, 


.43 


Q 


4C 31.63 


22 


03 


15 


.98 


+31 


45 


38. 


.3 


0. 


,295 


15. 


.58 


± 


0. 


05 f 


-25. 


,19 


L 


PG2214+139 


22 


17 


12 


.35 


+ 14 


14 


21. 


.2 


0, 


,066 


15, 


.05 


± 


0. 


.03' 


-22, 


,19 


Q 


PG2251+113 


22 


54 


10 


.50 


+ 11 


36 


39 


.0 





,326 


16 


.0 


± 





,l d 


-25. 


,02 


L 


PG2349-014 


23 


51 


56 


.13 


-01 


09 


13. 


.4 


0. 


,174 


15. 


.45 


± 


0. 


,20 j 


-24. 


,02 


L 



a Redshift given by NED 

h H a = 75, M = 0.3, n A = 0.7. 

C Q - radio quiet, L - radio loud: following Kellermann et al. (1994) where sources with a 
5GHz radio luminosity greater than 10 25 W Hz -1 (H = 50 km s _1 Mpc -1 , q — 0.5) were 
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referred to as radio loud and the weaker ones as radio quiet. 
d B magnitude from Giveon et al. (1999). 
e B magnitude from Surace & Sanders (2000). 
{ B magnitude from Neugebauer et al. (1979). 
K B magnitude from Schmidt & Green (1983). 
h B magnitude from Surace, Sanders, & Evans (2001). 
l B magnitude from McAlary et al. (1983). 
3 B magnitude = 15.45 (Simbad). 
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Table 2. Summary of Observations 



Name 


Obs. Date 
(UT) 


Exposure a 
(s) 


Filter 


FWHM 

(") 


Telescope 


Yt a a o c* i i oa 

PG0026+129 


OAA1 T A T 

2001 Jun 27 


2670 


TT 

ti 


0.13 


Gemini 


PG0050+124 


OAAA "\T — 1 T 

21)1)1) JNov 17 


/ 1 OAA\ 

(1200) 


JT 

ti 


A O O 

0.32 


Gemini 




OAAA 1\T 1 T 

2000 Nov 17 


(1620) 


K 


0.23 


Gcmmi 




OAAA AT— 1 A 

2000 Nov 19 


(1560) 


ti 


a i a 

0.16 


Gemini 




OAAA "\T 1 A 

2000 Nov 19 


(1200) 


K 


0.15 


Gcmmi 




OAAA 1\T OA 

2000 Nov 20 


/ O O A A \ 

(2300) 


H 


0.23 


Gemini 




OAAA "NT— OA 

21)1)1) JNov 2(J 


(1130) 


K 


A O 1 

I). 21 


Gemini 


PG0157+001 


OA A 1 O — 1 A 

2001 Sep 14 


a a rr a 

4050 


TT 


A 1 A 

0.19 


Gemini 


PG0804+761 


OAAA T~V Ol 

2000 Dec 21 


2910 


H 


0.15 


Gemini 




OAAA "P\ — — 01 

2000 Dec 21 


or'on 

2620 


K 


A 1 E 

0.15 


Gemini 


.rGl)ooo+77(J 


onno t?~v. 1 o 
2003 rcb 12 


1 a £?n 
1460 


TT 

ti 


n o o 
0.32 


Subaru 


I ■> An 1 J i O A A 

PG0844+349 


OAA1 ( \ j 1 A 

2001 Oct 10 


600 


TT 

ti 


0.38 


Gemini 


PG0923+201 


OA A 1 1 O 1 

2001 leb 21 


3930 


TT 

H 


0.20 


Gcmmi 




onm oo 
2001 reb 23 


/ E TAA\ 

(57UUJ 


J 


A 1 T 

U.17 


Gemini 


PG0953+414 


2000 Dec 20 


A CT A A 

4590 


K 


A 1 O 

0.13 


Gemini 




OAA1 T7 1 OO 

2001 leb 23 


(1380) 


J 

J 


0.23 


Gemini 


I ■> /"i 1 A A 1 i A f A 

PG1001+054 


OAAA T"\ o 

2000 Dec 26 


4140 


H 


0.20 


Gcmmi 


nrii nn^ i 1 OA 

PG1004+130 


OA AO T7 1 -. 1- 1 O 

21)1)3 rcb 12 


O A AA 

2400 


TT 

ti 


A 1 C 

0.16 


bubaru 


PG1116+215 


O A A O A O 1 

2002 Apr 21 


/ O A O A \ 

(3030) 


TT 

ti 


0.18 


Gcmmi 


Yt 1 1 1 A i 1 AA 

PG1119+120 


OAA1 T7 1 OO 

2001 teb 23 


/ O A /"* A \ 

(3960) 


TT 

ti 


0.15 


Gemini 


PG1126— 041 


OA A 1 T7* 1 1 

2001 leb 21 


4950 


TT 

ti 


0.15 


Gcmmi 


PG1211 + 143 


onno A 01 

2002 Apr 21 


/oata\ 

(2970) 


TT 

ti 


AIT 

0.17 


Gemini 


nr'ioifi i OA/1 

PG1229+204 


onno A oo 

2002 Apr 23 


(630) 


TT 

ti 


A O A 

I). 24 


Gemini 


PG1302 — 102 


2003 leb 12 


1 fOA 

1620 


TT 

ti 


Air 

0.15 


Subaru 


PG1307+085 


2003 Feb 12 


2420 


H 


0.15 


Subaru 


PG1309+355 


2002 Apr 23 


(2490) 


H 


0.22 


Gemini 


PG1351+640 


2002 Apr 21 


(2100) 


K' 


0.23 


Gemini 


PG1411+442 


2001 Feb 23 


(1800) 


H 


0.13 


Gemini 




2001 Feb 23 


(1620) 


J 


0.15 


Gemini 




2001 Feb 25 


3870 


H 


0.17 


Gemini 




2001 Apr 25 


2730 


K' 


0.13 


Gemini 


PG1426+015 


2001 Apr 24 


(4080) 


K' 


0.18 


Gemini 


PG1435-067 


2001 Feb 23 


(2100) 


J 


0.21 


Gemini 


PG1440+356 


2001 Jun 28 


(960) 


K' 


0.16 


Gemini 


PG1613+658 


2001 Feb 25 


4500 


K' 


0.30 


Gemini 




2001 Apr 25 


3420 


J 


0.18 


Gemini 
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Table 2 — Continued 



Name 


Obs. Date 


Exposure 3, 


Filter 


FWHM 


Telescope 




(UT) 


(s) 




(") 






2001 Apr 27 


1380 


H 


0.18 


Gemini 




2001 Apr 27 


1110 


K' 


0.28 


Gemini 


PG1617+175 


2002 Apr 21 


(4140) 


K' 


0.23 


Gemini 


PG1626+554 


2001 Apr 27 


1170 


K' 


0.18 


Gemini 




2001 Jun 27 


4530 


H 


0.16 


Gemini 


PG1700+518 


2001 Sep 12 


990 


H 


0.13 


Gemini 




2001 Sep 14 


3270 


H 


0.14 


Gemini 




2001 Oct 10 


1350 


K' 


0.15 


Gemini 


PG2130+099 


2001 Jun 28 


(3180) 


K' 


0.18 


Gemini 




2001 Oct 10 


1080 


H 


0.15 


Gemini 


4C 31.63 


2001 Jun 27 


3540 


H 


0.13 


Gemini 


PG2214+139 


2001 Sep 12 


2070 


H 


0.17 


Gemini 


PG2251+113 


2001 Sep 14 


2910 


H 


0.14 


Gemini 


PG2349-014 


2001 Sep 12 


3540 


H 


0.13 


Gemini 



a Values in parentheses indicate poor photometric conditions, hence the sen- 
sitivity may not scale as \ft 



Table 3. Near-Infrared Photometry 



Name 


Band 


2MASS a 


N79,N87 b < c Others 


This Work c Adopted Value d 






(mJy) 


(mJy) (mJy) 


(mJy) (mJy) 



PG0026+129 


H 


6.24 ± 0.2 


6.02 ± 1.0 




6.38 


6.38 


PG0050+124 


J 


16.69 ± 0.5 










PG0050+124 


H 


25 88 ± 1 










PG0050+124 


K' 


47.42 ± 1.0 










PG01 57+001 


H 


6 82 ± 2 


5 37 ±03 




6.16 


6.16 


PG0804+761 


H 


14.95 ± 0.2 


12 88 ± 6 




15.95 


15.95 


PG0804+761 


K' 


26.57 ± 0.2 


25.11 ± 2.0 




29.50 


29.50 


PG0838+770 


H 


3.11 ± 0.2 


3.09 ± 0.2 




2.94 


2.94 


PG0844+349 


H 


8.73 ± 0.1 


9.54 ± 0.5 




9.31 


9.31 


PG0923+201 


J 


3.12 ± 0.06 






> 3.07 


3.12 


PG0923+201 


H 


4.54 ± 0.2 






4.73 


4.73 


PG0953+414 


J 


3.40 ± 0.1 


3.39 ± 0.3 




> 2.73 


3.40 


PG0953+414 


K' 


6.52 ± 0.2 


7.76 ± 0.5 




6.72 


6.72 


PG1001+054 


H 


2.18 ± 0.2 


7.76 ± 0.5 


2.19 ± 0.1 c 


2.50 


2.50 


PG1 004+1 30 


H 


3.74 ± 0.3 


4.07 ± 0.4 


5.1 ± 0.2 f 


3.62 


3.62 


PG1116+215 


H 


8.84 ± 0.2 






> 10.54 


8.84 


PG1119+120 


H 


9.35 ± 0.1 


9.77 ± 0.3 




> 8.74 


9.35 


PG1126-041 


H 


15.78 ± 0.6 






15.49 


15.49 


_L VJT J_ Z( J_ J_ |^ _1_ 4: 


H 


Lzi.Ljkj —1— U . *J 


14.8 ± 1.0 




> 8 46 


12.23 


PG1229+204 


H 


8.67 ±0.1 


8.51 ±0.6 




> 6.23 


8.67 


PG1302-102 


H 


3.83 ±0.2 




3.74 ±0.2° 


3.99 


3.99 


PG1307+085 


H 


4.44 ±0.2 


4.32 ±0.4 




3.65 


3.65 


PG1309+355 


H 


4.06 ±0.1 






> 3.36 


4.06 


PG1351+640 


K' 


12.12 ±0.1 


10.96 ±0.8 




> 5.65 


12.12 


PG1411+442 


J 


8.03 ±0.1 


5.62 ±0.4 




> 7.19 


8.03 


PG1411+442 


H 


10.57 ±0.1 


8.32 ±0.8 




12.23 


12.23 


PG1411+442 


K' 


16.90 ±0.2 


17.38 ± 1.3 




16.8 


16.8 


PG1426+015 


K' 


24.59 ±0.4 


15.85 ± 1.1 




> 16.27 


24.59 


PG1435-067 


J 


2.76 ±0.1 






> 2.18 


2.76 


PG1440+356 


K' 


24.73 ±0.2 


25.12 ± 1.8 




> 17.56 


24.73 


PG1613+658 


J 


11.16 ± 0.2 


4.57 ±0.3 




7.87 


7.87 


PG1613+658 


H 


14.29 ±0.2 


6.31 ±0.6 




11.23 


11.23 


PG1613+658 


K' 


25.16 ±0.7 


10.72 ±0.8 




17.85 


17.85 


PG1617+175 


K' 


24.59 ±0.4 


12.30 ±0.9 




> 7.21 


24.59 


PG1626+554 


H 


3.70 ±0.3 






4.13 


4.13 
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Table 3 — Continued 



Name 


Band 


2MASS a 


N79,N87 b < c 


Others 


This Work 


Adopted Value d 






(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


PG1626+554 


K' 


6.15 ±0.1 






5.37 


5.37 


PG1700+518 


H 


6.44 ±0.2 


6.53 ±0.4 




6.19 


6.19 


PG1700+518 


K' 


10.69 ±0.1 


12.74 ±0.9 




9.38 


9.38 


PG2130+099 


H 


21.89 ±0.6 


15.37 ± 1.1 




18.74 


18.74 


PG2130+099 


K' 


36.64 ±0.5 


27.33 ±2.0 




> 19.88 


36.64 


4C 31.63 


H 


4.65 ±0.1 


6.17±0.5 




4.5 


4.5 


PG2214+139 


H 


13.82 ±0.2 




13.20 ± 0.8 s 


16.44 


16.44 


PG2251+113 


H 


4.38 ±0.2 


4.57±0.3 




4.48 


4.48 


PG2349-014 


H 


5.33 ±0.2 


5.75 ±0.3 


5.35±0.9 h 


5.58 


5.58 



a Flux enclosed within a 8" diameter disk. 

b Neugebauer et al. (1979), Neugebauer et al. (1987). Flux enclosed within a 5'/5 diameter disk 

c Flux enclosed within a 5'/5 diameter disk. Typical la uncertainty ±5%. Lower limits were 
derived in non-photometric conditions, and are accurate to ±5% (la level). 

d If available, the photometry done in this work is adopted. If not, 2MASS values are adopted. 

e Hyland & Allen (1982) 

f Sitko et al. (1982) 

sBalzano & Weedman (1981) 

h Rudy et al. (1982) 
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Table 4. Results of Host 2-D Fits 



Name/ 
Filter 



PG0026+129 H 
PG0157+001 H 
PG0804+761 H 
PG0804+761 K' 
PG0838+770 H 
PG0844+349 H 
PG0923+201 J 
PG0923+201 H 
PG0953+414 J 
PG0953+414 K' 
PG1001+054 H 
PG1004+130 H 
PG1116+215 H 
PG1119+120 H 
PG1126-041 H 
PG1211+143 H 
PG1229+204 H 
PG1302-102 H 
PG1307+085 H 
PG1309+355 H 
PG1351+640 K' 
PG1411+442 J 
PG1411+442 H 
PG1411+442 K' 
PG1426+015 K' 
PG1435-067 J 
PG1440+356 K' 
PG1613+658 H 
PG1613+658 J 
PG1613+658 K' 
PG1617+175 K' 
PG1626+554 H 
PG1626+554 K' 
PG1700+518 H 
PG1700+518 K' 
PG2130+099 H 



X r l/4 Xcxp 



1 


05 


1 


08 


3 


79 


4 


26 


1 


54 


1 


90 


1 


01 


1 


10 


4 


43 


5 


05 


2 


33 


3 


07 


1 


12 


1 


24 


1 


07 


1 


24 


o 


94 


o 


90 


o 


92 


o 


93 


o 


97 


1 


03 


1 


13 


1 


79 


1 


20 


1 


13 


11 


38 


10 


8G 


11 


61 


17 


00 


1 


28 


1 


28 


1 


97 


2 


46 


2 


20 


2 


46 


1 


42 


2 


15 


o 


98 


1 


28 


n 
u 


88 
oo 


n 
u 


88 


2 


08 


2 


05 


9 


20 


8 


77 


1 


48 


2 


26 


1 


53 


1 


52 





98 


1 


04 


2 


14 


1 


93 


4 


30 


G 


67 


5 


14 


5 


61 


G 


15 


13 


17 


1 


05 


1 


62 





89 





86 





93 





92 


5 


48 


5 


52 


2 


53 


2 


34 


4 


57 


3 


13 



1 


04 


[0.96/0 


04] 


3 


12 


[0.64/0 


36] 


1 


05 


[0.82/0 


18] 





95 


[0.68/0 


32] 


3 


58 


[0.78/0 


22] 


2 


31 


[0.69/0 


31] 


1 


15 


0.89/0 


11] 


1 


07 


[0.94/0 


06] 





88 


[0.32/0 

L / 


681 





91 


[0.98/0 


02] 





93 


0.48/0 


52] 


1 


06 


[0.65/0 


35] 


1 


13 


[0.00/1 


00] 


G 


57 


[0.70/0 


30 

J 


9 


66 


[0.64/0 


36] 


1 


30 


[0.46/0 


54] 


1 


41 


[0.65/0 


35] 


2 


00 


0.74/0 


26 

J 


1 


19 


[0.74/0 


26] 





97 


[0.92/0 


08] 





87 


[0.48/0 


52] 


2 


00 


[0.86/0 


14] 


8 


57 


[0.61/0 


39] 


1 


47 


[0.96/0 


04] 


1 


37 


[0.91/0 


09] 





96 


[0.92/0 


08] 


1 


65 


[0.84/0 


16] 


3 


95 


[0.90/0 


10] 


5 


15 


[0.99/0 


01] 


5 


97 


[0.98/0 


02] 


1 


06 


[0.98/0 


02] 





84 


[0.52/0 


48] 





92 


[0.79/0 


21] 


5 


38 


[0.80/0 


20] 


2 


34 


[0.00/1 


00] 


3 


28 


[0.00/1 


00] 
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Table 4 — Continued 



Name/ x r 2 i/4 xL P X r 2 i/4 +cxp a 



Filter 



PG2130+099 K' 5.25 

4C 31.63 H 1.07 

PG2214+139 H 1.79 

PG2251+113 H 1.04 

PG2349-014 if 4.31 



3.98 3.87 [0.42/0.58] 

1.49 1.06 [0.84/0.16] 

3.26 1.49 [0.76/0.24] 

1.02 0.95 [0.53/0.47] 

5.70 4.13 [0.82/0.18] 



a The numbers in brackets give the fraction of the 
host's luminosity originating from the bulge compo- 
nent (first number) and the disk component (second 
number) for the bulge+disk 2-D fit. 
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Table 5. Host Morphologies 3, 



Name 


T 

J 


•\^2 T~i b 

x -Aype 


A"' 


Visual Features 
Bar Arms Dist c 


/laoptcci lypc 


j-N OtCS 


PG0026+129 




( r l/4) 




no 


no 


no 


E 




PG0050+124 








no 


yes 


no 


p 




PG0157+001 




r 1 / 4 +exp. 




no 


yes 


str 


? 


host asymmetry 


PG0804+761 




r l/4 


r l/4 


yes 


no 


no 


P 


well- defined bar 


PG0838+770 




r l/4 




yes 


yes 


yes 


P 


twisted bar+ring, well-defined bar 


PG0844+349 




r l/4f 




no 


yes 


no 


Dp 




PG0923+201 


r l/4 


r l/4 




no 


no 


no 


E 




PG0953+414 


(cxp) 




? 


no 


no 


_ 




host too faint to be classified 


PG1001+054 




( r l/4) 




yes 


no 


no 


Dp 


well-defined bar 


PG1004+130 




r l/4 




no 


no 


no 


E 




PG1116+215 




exp 




no 


no 


yes 


Dp 


CC, host asymmetry 


PG1119+120 




r 1 ' 4 +exp. 




no 


no 


no 


B+D 


CC 


PG1126-041 




r 1/,4 +exp. 




no 


yes 


no 


B+D 


CC 


PG1211+143 




? 




no 


no 


_ 




host too faint to be classified 


PG1229+204 




r 1 / 4 +exp. 




yes 


yes 


no 


B+D 


well-defined bar 


PG1302-102 




r l/4 




no 


no 


no 


E 


CC 


PG1307+085 




r l/4 




no 


no 


no 


E 




PG1309+355 




r l/4 




no 


no 


no 


E 




PG1351+640 






? 


no 


no 


no 


? 




PG1411+442 


exp 


r 1 / 4 +exp 


r l/4 


no 


?g 


str 


? 


CC 


PG1426+015 






r l/4 


no 


?g 


str 


? 


CC 


PG1435-067 


r l/4 






no 


no 


no 


E 




PG1440+356 






exp 


no 


yes 


no 


Dp 




PG1613+658 


r l/4 


r l/4 


r l/4 


no 


?g 


str 


? 


CC 


PG1617+175 






r l/4f 


no 


no 


no 


E 




PG1626+554 




(cxp) 1 ' 


? f 


no 


no 


no 




bad PSF; unreliable 2-D fit 


PG1700+518 




r l/4 


exp 


no 


?g 


str 


? 


CC 


PG2130+099 




exp 


r 1//4 +cxp 


no 


yes 


no 


Dp 




4C 31.63 




r l/4 




no 


no 


no 


E 




PG2214+139 




r l/4 




no 


no 


yes 


E 


concentric shells 


PG2251+113 




? f 




no 


no 


no 




bad PSF; unreliable 2-D fit 


PG2349-014 




r l/4 




no 


?g 


str 


? 


CC 



a The radial profile in the central 1" to 2" radius of each host galaxy (200 pc to 2 kpc, depending on redshift) is usually not reliable, 
as PSF variations dominate this inner region. Classifications presented in this table therefore apply only to the region outside this 
central radius. 

b Bcst 2-D fit(s) correspond to the fit with minimum \ 2 m each band as listed in Table 6. Parentheses indicate 2-D fits with low 
confidence level (difference between X 2 i/4. an d Xexp ' s between 0.02 and 0.05). 

c Apparent degree of disturbance from visual inspection of the images: Disturbed ("yes"), Non disturbed ("no") or Strongly 
disturbed ("str"). 

d D p : disk present, E: Elliptical, B+D: Bulgc+disk. Type "disk present" indicates the presence of a disk component, but does not 
exclude a small bulge in the inner part of the host. The final adopted host morphology is based upon minimum \ 2 2-D fitting of 
images modified by visual inspection of the image for signs of features indicative of the presence of a disk component. 

e CC: Close (within lOkpc projected separation) compact companion (see text) 

f Bad PSF; 2-D fit is uncertain. 
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s Unclear whether features seen are arms or tidal debris. 
h No reliable 2-D fit available: does not appear in table 6. 
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Table 6. Host 2-D Fit Parameters for the Adopted Type. 



Name/ 
Filter 




Nuclear magn 
app. abs. 


Host 
app. 


magn 




m/" 2 


r e a 
(") 


e 




(rad) 


•r a 

m/" 2 


(") 




PA a 
(rad) 


PG0026+129 


H 


13 


05 


-25 


93 


14 


46 


-24. 


52 


22.16 


10.62 





21 





.01 


















PG0157+001 


II 


13 


41 


-25 


90 


13 


30 


-26. 


01 


18.79 


1.94 





.43 





60 


18 


.55 


3 


.75 





43 


1 


73 


PG0804+761 


H 


12 


.22 


-25. 


94 


14 


24 


-23 


92 














17 


.08 


1 


.74 





.28 





.82 


PG0804+761 


K' 


10 


67 


-27. 


49 


14 


44 


-23 


72 














18 


.02 


3 


.39 





.61 





.82 


PG0838+770 


II 


14 


31 


-24. 


48 


13 


.89 


-24 


90 














17 


26 


2 


.97 





59 


2 


.85 


PG0844+349 


II 


13 


35 


-23 


82 


13 


65 


-23 


52 














16 


.04 


1 


.38 





.25 





.53 


PG0923+201 


J 


14 


52 


-25 


16 


15 


46 


-24. 


22 


22.25 


5.73 





14 


1 


.18 


















PG0923+201 


II 


13 


53 


-26. 


15 


14 


.73 


-24 


95 


21.20 


4.53 





.07 





.79 


















PG0953+414 


J 


14 


39 


-25 


80 


15 


25 


-27 


77 


22.25 


3.58 





42 


3 


.07 


20 


95 


7 


.91 





40 


1 


.62 


PG0953+414 


K' 


12 


.42 


-24. 


94 


13 


86 


-26. 


33 


18.29 


0.95 





90 


1 


.55 


22 


02 


8 


.75 





.90 


1 


.98 


PG1001+054 


II 


14 


21 


-25 


06 


16 


68 


-22 


59 














19 


.59 


1 


61 





.12 





.56 


PG1004+130 


II 


13 


.84 


-26. 


41 


14 


12 


-26 


13 


22.00 


10.04 





.10 


2 


94 


















PG1116+21B 


II 


12 


.61 


-26. 


90 


15 


.96 


-23 


55 














17 


.50 





.99 





.33 





04 


PG1119+120 


H 


13 


.30 


-23 


32 


12 


.39 


-24 


23 


17.22 


0.91 





.43 


2 


.37 


16 


.07 


1 


.49 





.37 


1 


.38 


PG1126— 041 


H 


12 


.49 


-24. 


53 


12 


63 


— 24 


39 


17.94 


2.65 





.79 


1 


.21 


16 


.05 


1 


.65 





.52 


1 


00 


PG1211 + 143 


II 


12 


.37 


-25 


33 


14 


44 


-23 


26 


21.98 


4.99 





00 


2 


41 


20 


19 


7 


.49 





.51 





.85 


PG1229+204 


II 


13 


.56 


-23 


57 


12 


18 


-24 


95 


17.59 


1.57 





24 


1 


.51 


16 


.94 


3 


.38 





.59 


2 


67 


PG1302-102 


II 


13 


.74 


-26. 


88 


14 


.25 


-26. 


37 


19.89 


3.14 





46 





.81 


















PG1307+085 


II 


13 


97 


-25 


22 


14 


.71 


-24 


48 


21.38 


4.70 





.00 





64 


















PG1309+355 


II 


14 


.03 


-25 


57 


13 


71 


-25. 


89 


19.53 


3.15 





.18 


1 


.46 


















PG1351+640 


K' 


12 


.02 


-25 


84 


13 


77 


-24 


09 


20.83 


3.82 





.00 





64 


15 


.48 





.63 





.00 


3 


04 


PG1411+442 


J 


13 


48 


-24. 


43 


14 


.39 


-23 


52 


21.01 


6.32 





52 


1 


.28 


17 


20 





.82 





.55 


2 


.36 


PG1411+442 


H 


12 


.45 


-25 


46 


13 


.45 


-24 


46 


18.82 


1.30 





15 


1 


09 


17 


83 


2 


.85 





.56 


1 


.35 


PG1411+442 


K' 


11 


52 


-26. 


39 


13 


.54 


-24 


37 


19.84 


5.90 





.57 


1 


.28 


15 


99 





.62 





.83 


2 


50 


PG1426+015 


K' 


11 


.37 


-26. 


46 


12 


17 


-25 


66 


18.61 


5.33 





42 


2 


.45 


14 


.88 





.65 





.60 





36 


PG1435-067 


J 


14 


72 


-23. 


98 


15 


.25 


-23 


45 


21.13 


2.97 





00 


2 


.55 


















PG1440+356 


K' 


11 


51 


-26 


12 


12 


.80 


-24 


84 














14 


.63 


1 


.02 





18 


2 


.35 


PG1613+658 


J 


13 


56 


-25 


19 


13 


.76 


-24 


99 


20.98 


7.90 





22 





.82 


15 


99 





.30 





.90 


2 


14 


PG1613+658 


II 


12 


.90 


-25 


85 


12 


88 


-25 


87 


19.98 


6.94 





.25 





.95 


15 


.91 





.68 





46 


2 


91 


PG1613+658 


K' 


11 


.88 


-26 


99 


12 


11 


-26. 


64 


18.34 


4.00 





.14 





.88 


17 


.08 





.86 





.61 


2 


.77 


PG1617+175 


K' 


11 


12 


-27 


30 


13 


.40 


-25. 


02 


21.73 


14.34 





.00 





.08 


















PG1626+554 


II 


13 


.71 


-25 


12 


14 


33 


-25 


87 


20.34 


2.08 





.13 


1 


.69 


18 


53 


1 


.99 





.07 





.32 


PG1626+554 


K' 


12 


.96 


-24 


50 


13 


.11 


-25. 


72 


17.76 


0.30 





.01 


1 


.13 


15 


61 





.59 





.03 


1 


.30 


PG1700+518 


II 


13 


28 


-27. 


47 


14 


.62 


-26 


13 


20.85 


4.29 





.44 


1 


.98 


15 


.80 





.60 





.74 


1 


.59 


PG1700+518 


K' 


12 


04 


-28. 


71 


15 


.56 


-25. 


19 


inf 


2.11 





54 


1 


24 


15 


.87 


1 


.04 





81 


1 


70 


PG2130+099 


II 


12 


05 


-25 


08 


13 


.60 


-23 


53 














15 


.54 


1 


.25 





40 


2 


.27 


PG2130+099 


K' 


10 


76 


-26 


37 


12 


72 


-24 


41 














14 


42 


1 


.08 





.35 


2 


.24 


4C 31.63 


II 


13 


.68 


-27. 


09 


14 


25 


-26. 


52 


20.97 


6.47 





.37 





64 


















PG2214+139 


II 


12 


41 


-24. 


83 


12 


37 


-24 


87 


18.07 


2.72 





.07 


3 


.02 


















PG2251 + 113 


II 


13 


51 


-27 


51 


14 


.27 


-26. 


75 


19.77 


1.43 





34 





.21 


20 


.79 


1C 


i.OO 





.16 





56 


PG2349-014 


H 


13 


54 


-25 


93 


13 


57 


-25. 


90 


20.41 


6.19 





38 





22 


19 


32 


3 


;.60 





.51 


1 


67 



^Definitions of the parameters are given in equation 2. The bulge component is characterized by I e (effective surface brightness), r e 
(effective radius), e e (elongation) and PA e (position angle). The disk component is characterized by Ig (peak surface brightness), rg 
(effective radius), eg (elongation) and PAg (position angle). 

b Thc host type adopted in Table 5 docs not necessarily match the lowest x 2 value shown here, as morphological structures in the images 
is also taken into account. As discussed in the text, classification from \ 2 values alone is unreliable. In columns 2 (Nuclear Mad.) to 11 
(PAg), the results of the fit corresponding to the final host classification in Table 5 are given. For hosts that are unclassified (PG0953+414, 
PG1211 + 143, PG1411+442, PG1426+015, PG1613+658, PG1626+554, PG1700 + 518 and PG2349-014), the results of the bulgc+disk fits 
arc shown. 
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Table 7. Host Classification Statistics: iJ-band luminosities 21 





Number b 


Nucleus Mh 






Host M H 




£nucl./-^host 






med. 


ave. 


O'ave 


med. 


ave. 


&ave 


med. 


ave. 


Cave 


Elliptical 


10 


-26.04 


-25.93 


0.24 


-24.91 


-25.24 


0.27 


1.83 


2.21 


0.40 


Disk present 


7 


-25.08 


-25.20 


0.35 


-23.55 


-23.78 


0.26 


4.17 


6.58 


2.61 


Bulge+Disk 


3 


-23.57 


-23.81 


0.30 


-24.39 


-24.52 


0.18 


0.43 


0.62 


0.22 


Unknown 


7 


-25.85 


-25.99 


0.24 


-25.87 


-25.39 


0.28 


1.23 


2.16 


0.54 


Non-disturbed 


18 


-25.40 


-25.38 


0.25 


-24.47 


-24.62 


0.24 


1.93 


2.86 


0.56 


Disturbed 


3 


-24.83 


-25.40 


0.62 


-24.87 


-24.44 


0.36 


0.96 


7.84 


5.73 


Strongly disturbed 


6 


-25.88 


-26.02 


0.27 


-25.88 


-25.61 


0.26 


1.13 


1.68 


0.39 


Total 


27 


-25.57 


-25.52 


0.20 


-24.58 


-24.82 


0.19 


1.69 


3.15 


0.81 



a If no H band images were acquired, H magnitudes for the QSO and host galaxy were derived from J or if', 
assuming J — H = 0.94 and H — K = 0.97 for the QSO nucleus (average values derived from the photometry 
obtained in our sample). For the host galaxy, J — H = 0.72 — 0.4 x z and H — K = 0.22 + 2.0 x z were assumed, 
where z is the redshift. 

a Only objects for which a meaningful fit was obtained are included in this table: PG0050+124, PG0953+414, 
PG1211+143, PG1626+554 and PG2251+113 were excluded. This subsamplc of 27 objects is identical to the 
one shown in the top panel of Fig. 7 and in Fig. 8. 
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Table 8. Mid- and Far-Infrared Source Properties' 



Name 


f 12 i_Lm 


/25/J.m 


f 60 /j,m 




log Lir, 


Lm,/LBB h 




(mjy) 


(mjy) 


(mjy) 


(mjy) 


(Lq) 




PG0026+129 


22 ±4 


< 40 


< 27 


< 80 


10.93 ±0.04 


0.06 ±0.01 


PG0050+124 


549 ± 11 


1097 ± 20 


2293 ± 17 


2959 ± 51 


11.91 ±0.00 


1.71 ±0.15 


PG0157+001 


136 ± 28 


525 ± 55 


2376 ± 56 


2308 ± 127 


12.59 ±0.01 


1.46 ±0.13 


PG0804+761 


171 ± 22 


210 ± 25 


190 ± 34 


116 ±41 


11.62 ±0.02 


0.31 ±0.02 


PG0838+770 


34 ±5 


102 ± 7 


174 ±9 


159 ±21 


11.47 ±0.01 


1.79 ±0.09 


PG0844+349 


126 ± 29 


209 ± 35 


163 ± 41 


294 ± 97 


11.16 ±0.03 


0.47 ±0.04 


PG0923+201 


< 98 


150 ± 45 


340 ± 57 


1020 ± 184 


12.10 ±0.03 


0.64 ±0.05 


PG0953+414 


< 89 


< 107 


170 ± 26 


< 315 


11.80 ±0.06 


0.17 ±0.03 


PG1001+054 


34 ±6 


45 ± 11 


32 ± 9 


78 ±27 


11.39 ±0.03 


0.42 ± 0.03 


PG1004+130 


< 91 


< 150 


191 ± 42 


< 284 


11.88 ±0.06 


0.19 ±0.03 


PG1116+215 


152 ± 46 


187 ± 56 


< 219 


< 285 


12.06 ±0.04 


0.46 ± 0.06 


PG1119+120 


120 ± 40 


280 ± 53 


546 ± 51 


746 ± 126 


11.12 ±0.02 


0.70 ±0.05 


PG1126-041 


104 ± 19 


309 ± 35 


669 ± 26 


1172 ± 134 


11.33 ±0.01 


0.68 ±0.04 


Df 1 Q1 1 1 1 A Q 


ID l ± OU 


'XAI A- A A 
o4 ( it 44 


oZ { it OU 


4Uo it 14U 


1 1 kk _i_ n no 
11.00 it U.Uz 


n -U n no 
U.00 it V.uZ 


PG1229+204 


76 ±25 


246 ± 30 


183 ± 36 


317 ±95 


11.09 ±0.03 


0.60 ±0.04 


PG1302-102 


32 ±6 


< 180 


343 ± 68 


343 ± 68 


12.28 ±0.03 


0.25 ±0.03 


PG1307+085 


< 54 


< 153 


212 ± 42 


155 ±31 


11.46 ±0.05 


0.23 ±0.03 


PG1309+355 


60 ± 18 


102 ± 31 


< 140 


< 192 


11.76 ±0.04 


0.29 ±0.03 


PG1351+640 


172 ±5 


532 ± 6 


757 ±8 


1167 ±26 


11.81 ±0.00 


0.77 ±0.03 


PG1411+442 


115 ± 10 


160 ± 12 


162 ± 17 


140 ± 28 


11.38 ±0.01 


0.36 ±0.02 


PG1426+015 


127 ± 26 


195 ± 40 


323 ± 43 


346 ± 100 


11.48 ±0.02 


0.92 ±0.06 


PG1435-067 


50 ± 15 


< 126 


141 ± 52 


< 315 


11.34 ±0.04 


0.26 ±0.04 


PG1440+356 


101 ± 12 


208 ± 15 


651 ± 21 


1042 ± 62 


11.52 ±0.01 


0.73 ±0.07 


PG1613+658 


87 ± 12 


231 ± 14 


635 ± 19 


1090 ± 59 


11.96 ±0.01 


0.57 ±0.03 


PG1617+175 


65 ± 13 


71 ± 17 


102 ± 20 


< 252 


11.32 ±0.03 


0.31 ± 0.02 


PG1626+554 


33 ± 10 


39 ± 12 


< 70 


70 ±21 


11.16 ±0.04 


0.18 ±0.02 


PG1700+518 


111 ± 18 


222 ± 20 


466 ± 34 


441 ± 69 


12.68 ±0.01 


0.76 ± 0.04 


PG2130+099 


186 ±9 


380 ± 10 


479 ± 12 


485 ± 145 


11.39 ±0.01 


0.58 ±0.03 


4C 31.63 


54 ± 11 


106 ± 40 


193 ± 30 


212 ±82 


12.36 ±0.03 


0.42 ± 0.03 


PG2214+139 


61 ±7 


95 ± 12 


337 ± 11 


< 282 


11.00 ±0.01 


0.29 ±0.01 


PG2251+113 


< 36 


< 66 


< 67 


< 214 


< 12.04 


< 0.25 


PG2349-014 


< 119 


< 180 


271 ± 51 


290 ± 50 


11.75 ±0.05 


0.30 ±0.04 



a IRAS and ISO fluxes, in mjy, from Sanders et al. (1989), Haas et al. (2000), Polletta et al. (2000), and Haas et al. (2003) 
. When several independent measurements exist, we list the variance weighted mean value. Uncertainties arc 1 a; all upper 
limits arc 3 cr. 

b See text for definitions of Lir,, the "Infrared Bump" luminosity, and Lbb> the "Blue Bump" luminosity. 
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Table 9. Host Classification Statistics: FIR luminosities 





Number 


Median 3 


Median 3 


Median 3 


Mean [a] b 






M B (mag) 


log L m (L Q ) 


Lib/ Lbb 


Lib./ Lbb 


Elliptical 


10 


-23.88 


11.61 


0.27 


0.29 [0.14] 


Disk present 


8 


-22.52 


11.49 


0.53 


0.81 [0.56] 


Bulge+Disk 


3 


-21.63 


11.12 


0.68 


0.66 [0.04] 


Unknown c 


11 


-23.85 


11.77 d 


0.38 


0.54 [0.39] 


Non-disturbed 


21 


-23.16 


11.42 d 


0.42 


0.47 [0.35] 


Disturbed 


3 


-22.19 


11.47 


0.46 


0.85 [0.67] 


Strongly disturbed 


6 


-23.93 


11.85 


0.67 


0.73 [0.39] 


Unknown 


2 


-24.04 


11.67 


0.27 


0.27 [0.10] 


Total 


32 


-23.23 


11.52 d 


0.42 


0.54 [0.42] 



3 The median of an even number of values is the average of the two central values. 

b For the QSO with L m /L B B upper limit (PG2251+113), the 1 a L m /L B B upper limit 
is used to derive both the mean and its standard deviation. 

c Includes 6 hosts classified in Table 5 as "strongly disturbed" . 

d The median was determined without including PG 2251+113, for which there is only 
an upper limit on Lib.- 



